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The  purpose  of  this  thesis  was  to  demonstrate  the 
feasibility  of  the  moving-bank  multiple  model  adaptive 
estimation  algorithms  as  applied  to  flexible  spacestructure 
control.  Moving-bank  multiple  model  adaptive  estimation 
and  control  is  an  attempt  to  reduce  the  computational 
loading  associated  with  the  implementation  of  a  full-scale 
multiple  model  adaptive  estimator/controller.  The  results 
of  this  thesis  showed  that  the  performance  of  the 
moving-bank  controller  converged  to  the  performance  of  a 
benchmark  controller  (ideal  performance).  Also,  the 
moving-bank  controller  substantially  outperformed  a  fixed 
bank  controller  with  a  discretization  that  covered  the 
range  of  parameter  variation. 
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Abstract 

/ 

This  investigation  focused  on  the  use  of  moving-bank 
multiple  model  adaptive  estimation  and  control  (MMAE). 
Moving-bank  MMAE  reduces  the  computational  burden  of  MMAE 
by  implementing  only  a  subset  of  the  Kalman  filters  (9 
filters  versus  100  in  this  research)  that  are  necessary  to 
describe  the  system  to  be  estimated/controlled.  Important 
to  the  development  of  the  moving-bank  MMAE  are  the  decision 
logics  governing  the  selection  of  the  subset  of  filters. 

The  decision  logics  cover  three  situations:  initial 
acquisition  of  unknown  parameter  values;  tracking  unknown 
parameter  values;  and  reacqulsltlon  of  the  unknown  prar- 
meters  following  a  ^jump"  change  in  these  parameter  values. 

This  thesis  applied  moving-bank  MMAE  to  a  rotating  two 
bay  truss  model  of  a  flexible  spacestr ucture .  The  rotating 
two  bay  truss  approximated  a  space  structure  that  had  a  hub 
with  appendages  extending  from  the  structure.  The  mass  of 
the  hub  was  large  relative  to  the  mass  of  the  appendage. 

The  mathematical  model  was  developed  using  finite  element 
analysis,  transformed  into  modal  formulation,  and  reduced 
using  a  method  referred  to  as  singular  perturbations. 

X 


i  p^hj' jij 


w”.'  v-.  ■-.  ^n 


# 


Multiple  models  were  developed  by  assuming  that  variation 
occurred  in  the  mass  and  stiffness  of  the  structure. 

Results  showed  that  the  primary  cause  for  the  Increase 
in  the  performance  of  the  moving-bank  algorithm  over  past 
efforts  was  the  choice  of  the  measurement  covariance  noise 
matrix.  This  research  effort  also  showed  that  the  space 
discretization  level  can  play  an  Important  role  in  the 
performance  of  the  moving-bank  MMAC;  therefore,  care  must 
be  taken  in  determining  the  discretization  level.  Finally, 
this  research  effort  showed  that  the  moving-bank  MMAC 
provided  increased  control  performance  over  a  fixed-bank 
controller  with  a  discretization  that  covered  the  entire 
range  of  parameter  variation.  In  addition,  the  moving-bank 
controller  converged  to  the  performance  of  a  benchmark 
control ler . 
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MOVING-BANK  MULTIPLE  MODEL  ADAPTIVE  ESTIMATION 
APPLIED  TO  FLEXIBLE  SPACESTRUCTURE  CONTROL 


1 .  Introduction 


A  problem  exists  In  some  applications  of  estimation 
and  control  which  Is  caused  by  the  uncertalnity  of  para¬ 
meters  in  mathematical  models.  These  uncertain  parameters 
reduce  the  accuracy  which  can  be  expected  from  the  system 
model  and  the  algorithms  based  on  that  model.  Their  values 
can  be  constant  yet  unknown,  slowly  varying,  or  changing 
abruptly.  These  changes  in  parameters  often  necessitate 
the  Identification  of  parameters  within  the  mathematical 
model  and  the  changing  of  the  mathematical  model  during  a 
real-time  control  problem.  This  is  often  referred  to  as 
adaptive  estimation  and  control.  One  important  approach  to 
this  class  of  problems  is  multiple  model  adaptive  estima¬ 
tion  (MMAE)  and  control  (MMAC).  The  specific  models  in  the 
multiple  model  structure  can  either  be  static  ("fixed-bank” 
MMAE)  or  dynamically  redeclared  in  real  time  ("moving-bank" 
MMAE).  This  research  concentrates  on  two  areas.  First,  it 
explores  the  cause  of  the  moving-bank  MMAE  performance 
degradation  seen  in  Karnick's  191  research.  In  his 
research,  the  moving-bank's  estimate  of  the  true  parameter 
wandered  nonsystematlcally  throughout  parameter  space. 
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Second,  depending  on  the  cause  of  the  performance  degra¬ 
dation,  to  fine  tune  either  the  moving-bank  MMAE  or  the 
fixed-bank  MMAE  for  best  possible  performance  In  a  bending 
spacestructure  control  application. 


1.1.  Background 

In  a  system  adequately  represented  by  a  linear 
stochastic  state  model,  a  full  bank  MMAE  algorithm  can  be 
used  to  alleviate  the  problem  of  uncertain  parameters. 

This  method  consists  of  designing  a  Kalman  filter  for  each 
possible  parameter  value.  A  Kalman  filter  Is  a  recursive 
data  processing  algorithm  which  determines  the  optimal 
state  estimate  for  a  linear  system  with  dynamics  and 
measurement  noises  modeled  as  white  and  Gaussian.  It  is 
further  assumed  that  the  uncertain  parameters  can  take  on 
only  discrete  values;  either  this  is  reasonable  physically 
or  discrete  values  are  chosen  from  the  continuous  parameter 
variation  range.  The  MMAE  approach  results  in  a  bank  of  K 
Kalman  filters,  where  K  is  the  number  of  possible  parameter 
values  (see  Figure  1-1).  For  each  Kalman  filter,  a  condi¬ 
tional  probability  is  produced  that  assumes  its  parameter 
is  "correct*',  given  the  measurement  history.  These  condi¬ 
tional  probabilities  are  based  on  the  characteristics  of 
the  residuals  of  each  Kalman  filter  and  are  used  as  a 
weighting  factor  for  the  state  estimate  produced  by  each 


Kalman  filter.  The  Kalman  filter  based  on  the  best  esti¬ 
mate  of  the  "true"  parameter,  should  have  residuals  that 
most  nearly  match  the  anticipated  characteristics  of  the 
residuals  (as  the  result  of  the  filter-computed  residual 
covariance  matrix)  and  should  have  a  large  conditional 
probability,  in  contrast,  the  filters  that  correspond  to 
parameter  values  that  are  bad  estimates  of  the  "true" 
parameter  should  have  residuals  that  are  larger  than 
anticipated,  and  they  should  have  smaller  conditional 
probabilities.  The  adaptive  state  estimate  is  then  deter¬ 
mined  by  either  a  probabilistically  weighted  average  of  the 
outputs  of  each  Kalman  filter,  as  seen  in  Figure  1-1  (this 
is  called  the  Bayesian  form  of  the  MMAE  estimator),  or  by 
setting  the  overall  state  estimate  equal  to  the  state 
estimate  of  the  filter  with  the  highest  conditional 
probability  (called  the  maximum  a  posteriori,  or  MAP 
estimator)  120:365-3691. 

MMAE  has  been  successfully  implemented  in  several 
estimation  and  control  problems.  The  use  of  a  MMAE 
algorithm  to  track  airborne  targets  has  been  researched 
many  times.  Three  examples  of  this  research  are  given  in 
references:  ClOl,  C163,  and  C293.  The  MMAC  has  also  been 
used  in  controlling  fuel  tank  fires  C363,  addressing 
terrain  correlation  C30],  and  generating  estimators  for 
problems  in  which  large  initial  uncertainties  cause 
non-adaptlve  extended  Kalman  filters  to  diverge  1283. 


However,  one  basic  problem  of  the  full  bank  MMAE 
algorithm  is  the  number  of  Kalman  filters  which  must  be 
computed  simultaneously.  For  example,  if  a  system  had  two 
uncertain  parameters  and  if  each  of  these  parameters  can 
take  on  ten  different  values,  then  100  Kalman  filters  need 
to  be  Implemented.  As  can  well  be  Imagined,  this  is  an 
extremely  high  computational  burden  C19:1876D. 

Several  approaches  have  been  used  to  alleviate  the 
computational  burden  of  a  MMAE  C3:53.  One  method  uses 
Markov  processes  to  model  the  parameter  variation  Cl;25]. 

A  process  is  considered  Markov  if  its  present  parameter 
value  depends  only  on  the  single  previous  parameter  value 
Cl;4183  and  not  an  entire  history  of  values.  Other  methods 
include:  using  "pruning"  and/or  "merging"  of  "decision 
trees"  of  the  possible  parameter  time  history  C24;253, 
hierarchically  structuring  the  algorithms  to  reduce  the 
number  of  filters  C41,  and  a  method  in  which  the  filter  is 
initialized  with  a  coarse  parameter  space  discretization, 
but  after  the  filter  converges  to  the  "nearest"  parameter, 
the  filter  is  rediscretized  using  a  simple  direct 
method  C153. 

To  reduce  this  burden  for  the  example  of  2  parameters 
with  10  possible  values  for  each,  a  moving  bank  of  Kalman 
filters  consisting  of  only  nine  filters  can  be  developed 
(see  Figure  1-2).  These  Kalman  filters  should  be 


PARAMETER  a  ^ 

Q  used  Kalman  filter 
•  unused  Kalman  filter 

X  current  best  estimate  of  the 
true  parameter  value 
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associated  with  the  parameter  values  that  most  closely 
surround  the  estimated  "true"  parameter  value.  When  the 
"true"  parameter  value  lies  within  the  region  of  the  bank 
of  Kalman  filters,  the  moving-bank  MMAE  algorithm  behaves 
essentially  the  same  as  the  full  bank  MMAE  algorithm. 
However  when  the  "true"  parameter  value  lies  outside  of  the 
bank  or  even  inside  the  bank  but  near  its  boundary,  the 
bank  must  be  moved  or  expanded  (see  Figure  1-3).  A  bank  is 
moved  or  expanded  by  changing  the  parameter  points  to  which 
the  nine  Kalman  filters  in  the  bank  correspond  tl9:18761. 

In  the  past  three  years,  research  has  been  conducted 
in  the  area  of  moving-bank  MMAE  algorithms  to  reduce  the 
amount  of  computational  loading  needed  to  estimate  both  the 
state  variables  and  uncertain  parameters  associated  with 
the  system  model.  Initial  research  conducted  by  Hentz  and 
Maybeck  [6;193  showed  that,  for  a  simple  but  physically 
motivated  system,  the  moving-bank  MMAE  algorithm  performed 
as  well  as  the  full-bank  MMAE  algorithm,  but  with  an  order 
of  magnitude  less  computational  loading.  They  also  showed 
that  the  system  performance  could  be  enhanced  by  starting 
the  bank  with  a  large  discretization  to  allow  for  Initial 
parameter  acguistlon  and  reducing  the  discretization  as  the 
accuracy  of  the  parameter  estimate  Increases.  In  the 
follow-on  research  effort  conducted  by  Flllos  C33,  he 
showed  that  the  uncertain  parameters  associated  with  the 


system  model  can  be  determined  using  ambiguity  function 
analysis.  The  research  conducted  last  year  by  Karnlck  and 
Maybeck  showed  that  a  large  discretized  fixed-bank  MMAE 
algorithm  performed  at  least  as  well  as  the  moving-bank 
MMAE  algorithm  for  a  more  complex  system  C93.  Furthermore, 
the  moving  bank  never  converged  to  a  consistent  parameter 
value.  Instead  the  moving  bank  wandered  nonsystemat Ically 
throughout  the  parameter  space.  The  fixed-bank  MMAE  algo¬ 
rithm  included  nine  Kalman  filters  similar  to  the  moving- 
bank  MMAE  algorithm;  however,  it  used  coarser  discretiza¬ 
tion  and  was  not  allowed  to  move  t9D. 

1.2.  Problem 

The  full  bank  MMAE  algorithm  is  computationally  too 
burdensome  for  most  applications.  The  moving-bank  MMAE  was 
shown  to  be  successful  for  a  physically  motivated  but 
simple  system  C6].  However,  the  moving-bank  MMAE  algorithm 
was  shown  to  have  not  performed  as  well  as  a  fixed-bank 
MMAE  algorithm  for  a  more  complex  space  structure  appli¬ 
cation  that  required  adaptive  estimation  and  control  C93. 

1.3.  Scope 


This  research  applies  to  problems  associated  with 
large  space  structures.  The  large  space  structure  model 
used  In  this  research  was  a  two-bay  truss  (developed  by 


Drew  Karnlck  in  his  thesis  C93)  attached  to  a  hub  (see 
Figure  3-2)  which  allowed  motion  in  two  degrees  of  freedom 
(x-y  plane),  but  did  not  allow  translational  motion  of  the 
hub.  The  truss  was  100  inches  long  and  18  Inches  high.  In 
addition,  non-structural  masses  were  added  to  the  truss 
model  for  two  purposes.  First,  they  allow  time  varying 
mass  problems  to  be  modeled,  such  as  the  mass  associated 
with  fuel  tanks  being  depleted.  Second,  if  the  nonstruc- 
tural  masses  were  large,  relative  to  the  structure  mass, 
they  produce  the  low  frequency  structural  model  associated 
with  large  space  structures.  The  space  structure  model  was 
described  in  terms  of  mass  and  stiffness  matrices  obtained 
from  a  finite  element  analysis. 

Two  uncertain  parameters  were  investigated,  the 
non-structural  mass  and  the  stiffness  matrix  within  a 
linear  system  model  structure.  These  parameters  were 
chosen  because  they  physically  appear  in  the  system  dynam¬ 
ics  equation.  However,  "nonphysical"  parameters  could  be 
used  if  found  to  produce  better  estimates.  The  dynamics 
and  measurement  noise  characteristics  were  assumed  known 
and  modeled  as  white  Gaussian  processes. 

1.4.  Approach 

This  research  concentrated  on  two  areas.  The  first 


was  to  determine  the  cause  of  the  moving-bank  MMAE 


performance  degradation  seen  in  Karnick's  research  and  the 
second  was  to  fine  tune  either  the  moving-bank  MMAE  or  the 
fixed-bank  MMAE,  depending  on  the  cause  of  the  performance 
degradation.  In  previous  thesis  efforts,  Hentz  C6I,  Filios 
C33,  and  Karnick  C9I  developed  extensive  software  programs 
which  were  also  used  in  this  research.  This  software 
Included  (1)  a  preprocessor,  which  created  a  parameter 
space  that  was  utilized  in  a  Monte  Carlo  simulation  and 
ambiguity  function  analysis  (to  be  discussed  at  length 
subsequently),  (2)  a  processor,  which  simulated  the  moving- 
bank  or  fixed-bank  multiple  model  adaptive  estimator, 
performed  Monte  Carlo  simulation  runs  and  generated  data 
for  each  run,  (3)  a  postprocessor,  which  computed  the  means 
and  variances  of  variables  of  interest  and  then  generated 
the  plots  of  statistics  for  the  Monte  Carlo  simulation,  and 
(4)  a  program  that  computed  the  ambiguity  functions  and 
generated  their  plots. 

Karnick  showed  that,  although  the  use  of  a  moving 
bank  may  provide  enhanced  state  estimation  performance, 
similar  performance  could  have  been  obtained  from  a  fixed 
bank  estimator  with  a  coarse  space  discretization  that 
covered  the  range  of  parameter  variation.  The  problem  that 
the  moving  bank  had  with  parameter  estimation  could  have 
been  caused  by  the  dynamic  driving  noise  (Q)  and  the  mea¬ 


surement  noise  (R)  values,  system  complexity,  or  unmodelled 
effects . 


After  the  appropriate  Q  and  R  values  were  determined, 
then  the  first  bending  mode  was  extracted  from  the  system 
model  and  the  truth  model.  This  was  done  so  that  an 
attempt  could  be  made  to  duplicate  Hentz's  work  C6l  for  the 
large  space  structure  system  model.  Hentz  showed  that,  for 
a  two-state  (i.e.  one  bending  mode)  filter  model  matched 
against  a  two-state  truth  model,  the  center  of  the  moving 
bank  would  move  in  the  direction  of  the  true  parameter. 

The  moving-bank  MMAE  algorithm  did  duplicate  Hentz's  work 
and  that  performance  did  not  degrade  substantially  when  the 
moving-bank  MMAE  algorithm  was  applied  to  the  entire  space 
structure  model  (where  the  space  structure  model  and  the 
truth  both  have  a  rigid  mode  and  two  bending  modes).  This 
suggests  that  the  benefits  gained  with  the  moving-bank  MMAE 
algorithm  do  not  degrade  as  the  system  complexity 
increases . 

Finally,  the  issue  of  robustness  in  the  face  of 
unmodelled  effects  could  be  addressed  by  augmenting  states 
to  the  truth  model  and  then  reevaluating  the  moving-bank 
MMAE,  based  on  the  three-bend Ing-mode  space  structure 
model,  against  the  new  truth  model.  By  examining  the 
eigenvalues  and  frequencies  of  the  space  structure  model 
before  state  reduction  was  performed  (see  Table  3-2),  a 
truth  model  with  16  states  would  be  recommended.  However, 
this  issue  was  not  addressed  in  this  effort.  Since  the 
problem  with  the  moving  bank  was  due  to  large  Q  and  R 


values,  the  decision  was  made  to  conduct  the  remaining 
research  using  a  moving-bank  MMAE  with  sufficiently  small  Q 
and  R  values. 


Following  this  decision,  the  space  discretization 
which  produced  the  best  state  estimate  was  determined.  An 
appropriate  space  discretization  level  was  needed  so  that 
the  moving-bank  algorithm  could  distinguish  between  incre¬ 
mental  steps  in  parameter  space.  In  previous  research 
efforts,  the  space  discretization  was  chosen  by  intuition. 
However,  investigation  was  needed  into  the  effect  of  var¬ 
ious  discretizations  of  the  parameter  space.  Insight  into 
the  parameter  space  discretization  could  be  acquired  from 
three  sources.  The  first  source  of  insight  could  be  ambi¬ 
guity  function  plots  [22;3;93.  Slowly  varying  ambiguity 
function  plots  versus  parameter  variations,  l.e.  flat, 
illustrate  that  the  parameter  value  associated  with  the 
plots  was  difficult  to  estimate,  but  they  also  Illustrate 
that  the  parameter  has  a  small  effect  on  the  state  estimate 
and  therefore  did  not  need  to  be  estimated.  In  contrast, 
quickly  varying  ambiguity  function  plots,  i.e,,  containing 
a  sharp  peak  and  valleys.  Illustrate  that  the  parameter 
could  be  estimated  and  does  substantially  affect  the  state 
estimate.  The  faster  the  ambiguity  function  plots  varied, 
the  tighter  the  space  discretization  needs  to  be.  However, 
Karnick  C93  experienced  extreme  difficulties  in  analysing 
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the  ambiguity  function  plots  because  they  consisted  of 
multiple  peaks.  The  multiple  peaks  made  determining  the 
size  of  the  space  discretization  level  difficult  if  not 
impossible.  The  second  source  of  Insight  could  be  the 
values  of  the  dynamics  driving  noise  strength  (Q)  and  mea¬ 
surement  noise  covariance  (R)  C243.  Large  Q  and  R  values 
indicated  a  need  for  large  space  discretization  because 
such  large  Q  and  R  values  would  mask  differences  in  the 
residuals  of  filters  in  a  tightly  spaced  discretization. 
This  difference  in  residuals  was  critical  for  the  operation 
of  a  MMAE .  The  third  source  of  insight  could  be  a  sensi¬ 
tivity  study  conducted  by  giving  the  filter  an  Incorrect 
parameter  value  and  then  observing  the  state  estimate  C91. 
The  more  inaccurate  the  state  estimate  becomes,  for  a  given 
size  parameter  variation,  the  tighter  the  space  discreti¬ 
zation  needs  to  be  for  that  parameter.  This  technique  was 
used  because  of  the  multiple  peak  problem  experienced  with 
the  ambiguity  function  plots. 

Following  the  space  discretization  study,  a  sensi¬ 
tivity  study  was  performed  to  determine  the  upper  bounds  on 
the  dynamics  driving  noise  (Q)  and  measurement  covariance 
noise  (R).  These  bounds  determine  the  precision  the  Q  and 
R  values  must  have  in  order  for  the  moving-bank  algorithm 
to  perform  correctly.  This  study  differs  from  the  previous 
Q  and  R  sensitivity  study  by  determining  the  actual  bound 
on  the  precision  of  the  Q  and  R  values,  while  the  previous 
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study  only  wanted  to  choose  Q  and  R  values  lower  than  the 
bound  values.  In  conducting  this  study,  the  R  values  were 
held  constant  while  the  Q  values  were  increased  one  at  a 
time  by  orders  of  magnitude  until  the  adaptive  mechaniza¬ 
tion  began  to  degrade.  Then  the  Q  values  were  held  con¬ 
stant  while  the  R  values  were  varied  one  at  a  time  by 
orders  of  magnitudes  until  the  adaptive  mechanization  began 
to  degrade. 

Following  the  Q  and  R  sensitivity  study,  a  study  was 
conducted  on  the  controller  performance  sensitivity  to 
altered  cost  function  definition.  In  the  research  Karnick 

performed,  the  quadratic  cost  function  used  to  determine 

* 

the  full-state  feedback  gain  (G  )  had  no  cross  terms. 

— c 

However,  if  control  influence  is  desired  on  the  states 
over  the  entire  sample  period  then  cross  terms  can  be 
non-zero.  Therefore,  in  this  study,  the  cross  weighting 
matrix  (S)  was  evaluated  C23:763  and  compared  to  the  state 
weighting  matrix  (X)  and  the  control  weighting  matrix  (U). 
The  comparison  showed  that  the  appropriately  chosen  cross 
weighting  matrix  was  small  compared  to  the  X  and  U  matri¬ 
ces.  This  was  mostly  due  to  the  very  small  sample  period 
time  as  compared  to  the  system  characteristic  times. 
Therefore,  the  cross  weighting  matrix  was  not  added  to  the 
quadratic  cost  function. 

Following  this  decision,  the  X  and  U  matrices  were 
tuned  to  provide  for  small  rms  values  on  the  true  states 


while  not  driving  the  control  Inputs  beyond  their  maximium 
value.  The  initial  X  and  U  matrices  were  equal  to  the  X 
and  U  matrices  used  by  Karnlck.  The  tuning  process  was 
accomplished  by  holding  the  U  values  constant  while  varing 
the  X  values  one  at  a  time,  then  the  X  values  were  held 
constant  while  varing  the  U  values  one  at  a  time. 

After  the  X  and  U  values  were  tuned,  a  study  was 
conducted  on  transitioning  from  nonadaptlve  controller 
gains  to  adaptive  controller  gains  based  on  parameter 
identification.  This  study  was  conducted  because  if 
adaptive  control  inputs  were  supplied  to  the  system  before 
an  accurate  estimate  of  the  true  parameter  was  obtained, 
the  control  Inputs  could  drive  the  system  unstable. 
Therefore,  a  nonadaptlve  controller  was  designed  that  was 
very  robust  to  parameter  variation  at  the  expense  of  con¬ 
troller  performance  so  that  the  controller  inputs  would  at 
least  maintain  stability  regardless  of  the  estimate  of  the 
true  parameter.  However,  the  moving  bank's  estimate  of  the 
true  parameter  was  locking  onto  the  actual  true  parameter 
so  quickly  that  the  nonadaptlve  controller  was  not 
required.  Therefore,  there  was  no  transition  from  non- 
adaptive  control  to  adaptive  control  in  this  research. 

The  final  study  was  a  comparison  between  three 
controllers  (a  moving-bank  controller,  a  fixed-bank  con¬ 
troller,  and  a  benchmark  controller).  Ideal  performance 


was  created  using  a  benchmark  controller.  In  the  benchmark 
controller  only  one  filter  was  matched  against  the  truth 
model  and  this  filter  was  given  artlfical  knowledge  of  the 
true  parameter.  The  fixed-bank  controller  was  designed 
identically  to  the  moving-bank  controller  except  the  nine 
implemented  filters  were  not  allowed  to  move  in  parameter 
space.  The  performance  criteria  used  in  the  controller 
study  was  the  mean  and  standard  deviation  on  the  true 
states . 

1 . 5  Summary 

The  remaining  chapters  in  this  thesis  cover  the 
following  areas.  Chapter  II  develops  the  algorithms  for 
the  moving-bank  MMAE  and  MMAC .  Chapter  III  develops  the 
large  spacestructure  model.  Chapter  IV  describes  the 
simulations  used  in  this  thesis  and  gives  a  brief 
description  of  how  the  simulations  are  to  be  used  in  the 
research.  Chapter  V  details  the  results  of  the  research. 
Finally,  Chapter  VI  contains  the  conclusions  and 


recommendations . 


2.1.  Introduction 


This  chapter  develops  the  algorithms  for  the 
fullscale  and  moving-bank  Bayesian  Multiple  Model  Adaptive 
Estimator.  First,  the  full-scale  model  Is  developed.  This 
is  then  modified  for  the  moving  bank  case.  Finally,  the 
Ambiguity  Function  analysis  is  developed. 

2.2.  Bayesian  Multiple  Mode  1  Adapt ive  Est imat ion 
Algor ithm  Development 

Development  of  the  full-scale  Bayesian  Multiple  Model 
Adaptive  Estimation  (MMAE)  algorithms  is  presented  in  this 
section.  For  a  more  rigorous  development,  the  reader  is 
directed  to  reference  122:129-1363. 

Let  the  system  under  consideration  be  discrete  and 
described  by  C3;  6;  21;  93: 

x(t,^,)  =  *(t, ,t, )x(t, )  +  B.(t.)u(t,)  +  G.(t,)w  (t.) 

z(t, )  =  H(t, )x(t, )  +  v(t,  )  (2-1) 

^  i  1  sc  1  »=  1 

where  denotes  a  vector  stochastic  process  and: 

x(t,)  =  n-d Imens i onal  state  vector 
^  1 


':yrrr»^Tr-frif:v-^-<mvTr^r.'\r_ 
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u(t^  ) 


=  state  transition  matrix 
=  r-dlmensional  known  input  vector 
=  control  input  matrix 


W(3(ti)  =  s-dlmenslonal  white  Gaussian  dynamics 
noise  vector 


noise  input  matrix 


z(t.)  =  m-dlmensional  measurement  vector 

JC  1 

H(tj^)  =  measurement  matrix 

v(ti)  =  m-dlmenslonal  white  Gaussian 
measurement  noise  vector 


and  the  following  statistics  apply: 


ECWj(tj>J=0 

ECv(t, ) ]  =  0 

JC  i  — 

ECv(t,  )v'^(t  .  )1  =  R(t,  )j,  ,  , 

«  1  ~  j  1  1 J 

where  6  is  the  Kronecker  delta  function.  It  is  also 
assumed  that  x(t-),  w,  ( t . ) ,  and  v ( t , )  are  Independent  for 
all  times  t^^. 

All  of  the  following  matrices  within  Equation  (2-1) 
are  functions  of  system  parameters:  ♦,  G^,  Q^,  H,  and 

R.  Now  suppose  these  system  parameters  have  values  that 
are  not  known  with  absolute  certainty  or  that  can  vary.  To 
address  this  problem,  let  a  be  the  uncertain  p-dlmensional 
parameter  vector  which  is  an  element  of  A,  where  A  is  a 
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p 

subset  of  R  (in  this  case  R  stands  for  the  set  of  real 

p 

numbers  and  R  Is  real  Euclidean  p-dimensional  space). 

This  parameter  vector  may  be  uncertain  but  constant.  It  0iay 
be  slowly  varying,  or  It  may  undergo  jump  changes.  The 
purpose  of  a  Bayesian  estimator  conceptually  Is  to  compute 
the  following  conditional  density  function  [223: 

),a|Z(ti  li)  *  fx(ti  )|  a,Z(ti  * 

^aiZ(ti) (alii  )  (2-2) 

where  Is  the  vector  of  measurements  from  t^  to  t^. 


Z(t,)  *  Cz’'(t,  ),z‘^(t,  ,  .  .,z'^(tn)3'^ 

The  second  term  on  the  right  side  of  Equation  (2-2) 
can  be  further  evaluated  r223: 

^a  |Z(ti )  (§.lli  )  =  ^|.||(ti  ),Z(ti-i)<5.l2.i'Ii-l) 

^a,|.(  1 1 )  I  Z.(  1 1-1 )  (3.^11 1  li  -1 1 

^j^(ti)  l(ti_i)(2.i  li-l) 


^j|(ti  )|  a,l(ti_i )  (Z.l|a,Ii-i)fa|Z(ti_i)  (a|Zi-i) 

r^t(ti  )|  a,Z.(ti_i)  (2lla,Zi-i)£g^l2(ti-i)  (a.lli-i)da. 


Conceptually,  Equation  (2-3)  can  be  solved 


recursively,  starting  from  an  a  priori  probability  density 
function  of  f^(al,  olncc  f  ,|  j,.  13 

Gaussian  with  a  mean  of  H(t^)^(t^)  and  covariance 
CH(  tj  )P{  t‘ )h'^(  t^  )+R(  t^  )  1 ,  where  t(t^)  and  P(t")  are  the 
conditional  mean  and  covariance  respectively  of  x(t,)  just 
prior  to  the  measurement  at  t^,  assuming  a  particular 
realization  a  of  a  C220. 

Using  the  conditional  mean,  the  estimate  of  x(t^) 
becomes  [ 22 J : 


ECx(ti ) I Z(ti  ) 


~  ^1^  x  { t  i  )  I  Z  ( t  i ^  ^  l^i  * 


fx(ti),a  (Z(ti)  >^3dx 

R  ^  (2-4) 


Ei:x(ti )  I  z(ti )  =  ^x(ti)  ||.,Z(ti)  ) 

^a  |Z(ti)  <§.  li 

-i  -  ^x(ti)|a,Z(ti)(2iia/Zi)^3fg|Z(tj)(§IZi)^ 

(2-5) 


where  the  term  in  brac)?ets  is  the  estimate  of  x(t,)  based 
on  a  particular  value  of  the  parameter  vector.  This  Is 
the  output  of  the  Kalman  filter  based  on  that  realization 
of  the  parameter  vector.  When  a  is  continuous  over  A,  this 
requires  an  infinite  number  of  filters  in  the  banit .  To 
reduce  the  number  of  filters,  the  parameter  space  is 


usually  discretized,  yielding  a  finite  number  of  filters. 
The  Integrals  over  A  In  Equations  (2-4)  and  (2-5)  then 
become  summations.  Defining  probability 

that  the  elemental  filter  Is  correct,  conditioned 

on  the  measurement  history,  it  can  be  shown  by  a  method 
analogous  to  the  development  for  Equation  (2-3)  that 
satlsf ies  C  221 : 


^ j|.(  t  i  )  1 1.,  |.(  t  i  _  1 )  <  li  I  2.1  - 1 )  * P)i  <  ^  1  - 1 ) 

Pk{ti)  =  - -  -  (2-6) 

K 

j?i  ^|.(ti)  |a,l(ti_i)  <£l  ia-j/2.1-1  )*Pj  (ti-i ) 


K 

x(tt)  =  ECx(ti  )( Z(ti  )  *  iil  =  53  >  ‘2-7) 

j  =  l 

where  a  C  C  a  ^  ,  a2  ,  .  .  . ,  aj^  l  and^(tj)  is  the  mean  of  2c(t^) 
conditioned  on  a  =  and  i.(t^)  *  l.e.  the  output  of 

the  k^^  Kalman  filter  In  the  bank,  based  on  the  assumption 
Plctorlally,  the  algorithm  appears  as  In 
Figure  1-1. 

The  probability  weighting  factors  for  each  Kalman 
filter  are  calculated  from  Equation  (2-6)  1223,  where 

^i^(ti  )  |fi.,|,(ti-i )  ‘ii  lS.kf2.i-l) 


1/1  (2n)"^2  )|  1/2  {  expC-(l/2)Ex(ti  )Axl(ti  )r)c(ti  )3 


Ak<ti)  =  Hk(ti)P,c(tI)H;^(ti)  +  R)((ti) 
LK(ti)  =  Zi  -  Hn(ti)i,j(tl) 

m  =  number  of  measurements 


Both  the  residual  covariance  Aj^(t^)  and  the  residual  r|^(t^) 
itself  are  readily  available  from  the  elemental  filter. 

The  estimate  of  the  parameter  and  the  covariance  of  the 
parameter  are  given  by  C223: 


a(tp  =  ECai  Z(ti  )  = 


Zil  a  falZ(ti)  (§l?i 

-/-OO  ^ 

=/  at  S  Pk(ti  )i(a  -  ak)3^ 

J-OCI 


(2-9) 


E’  Ca 


a(ti  )  iCa 
K 

=  k?i^^ 


-  a(ti )lT|^(ti )  =  Zil( 

-  a(ti)JCaj{  -  'a(ti)]’r  Pk(ti) 


(2-10) 


The  covariance  of  the  state  estimate  Is  given  by  C22]: 


P(tt)  =  ElCx(ti)  -  x(tt) JC^(ti )  -  x(tt)lT( z(ti  )  =  Zij 
=  f  lx  -  x(tt)lCx  -  i(tf)DTf^j,.^  J  ,2(ti)  (xlZ-i 

y.OO  ^ 

K  ^op 

=  k  =  l  Pk^^l^  /  ~  ^(ti)!"^ 

«/-00 

( 1 1 )  1  z  ( 1 1 )  ( ?^  1  iJt  f  ^1 )  111 


=  )C  =  1  Pk(ti)CPk(tt)  +  C^(tt)  -  x(tt)D* 

-  x;(t|-)3T3  (2-11) 

where  is  the  covariance  of  the  state  estimate  of  the 

elemental  filter. 

The  Bayesian  MMAE  has  been  shown  to  be  optimal  and  to 
converge  if  the  true  value  of  the  parameter  Is  nonvarying 
[51.  Convergence  for  this  case  occurs  when  the  probability 
associated  with  one  elemental  filter  Is  essentially  one  and 
the  probability  associated  with  all  other  elemental  filters 
Is  essentially  zero.  The  MMAE  converges  to  the  elemental 
filter  with  parameter  value  equal  to,  or  most  closely 
representing,  the  true  parameter  set,  as  defined  in  [53. 

There  are  no  theoretical  results  available  for 
varying  parameters  C3:18;  6:8;  9:203.  The  fact  that  the 
filter  can  converge  to  one  filter  for  a  non-varying  true 
parameter  value,  does  give  reason  for  some  concern.  For 
example.  If  the  true  parameter  value  is  varying  rather 
slowly,  the  algorithm  may  assume  one  filter  is  correct  with 
probability  essentially  equal  to  one.  However,  the  true 
parameter  value  may  eventually  become  significantly  dif¬ 
ferent  from  the  value  Initially  estimated  by  the  filter 
C6:93,  resulting  in  filter  divergence. 

Another  possibility  Is  that  the  algorithm  may  con¬ 


verge  and  lock  onto  the  "wrong"  filter.  The  filter  Is,  to 


probabilities  are  rescaled  so  that  their  sum  would  remain 
equal  to  unity. 


2.3,  Moving-Bank  Algor ithm  Development  C3:22-333 

The  MMAE  generates  a  computational  burden  that  is  too 
large  for  most  practical  applications  C3;  7;  19;  93. 

Maybeck  and  Hentz  demonstrated  that  the  full  bank  of 
filters  could  be  replaced  by  a  subset  of  filters  based  on 
discrete  parameter  values  "closest"  to  the  current  estimate 
of  the  parameter  vector.  The  probability  associated  with 
non- implemented  filters  is  conceptually  set  to  zero  while 
the  probability  weightings  are  distributed  among  the 
implemented  filters.  As  the  parameter  set  estimate 
changes,  filters  that  are  "closer"  to  the  new  parameter 
estimate  are  implemented  while  those  "farthest"  away  are 
removed.  Maybeck  and  Hentz  also  investigated  changing  the 
discretization  levels  of  the  moving  bank  model.  During  the 
acquisition  stage,  the  Implemented  filters  are  set  to  a 
coarse  discretization,  then  changed  to  finer  discretiza¬ 
tions  as  the  parameter  estimate  Improves.  Therefore,  the 
implemented  filters  do  not  necessarily  occupy  adjacent 
discrete  points  in  the  parameter  space,  as  would  be  used  in 


the  full  bank  MMAE. 


2.3.1. 


Average  C3;  6;  93.  The  outputs  of 


each  elemental  filter  of  the  moving  bank  estimator  are 
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weighted  and  summed  In  the  same  manner  as  Equations  (2-6) 
and  (2-7);  however,  only  the  Implemented  filters  In  the 
moving  ban)t  are  summed.  If  J  filters  are  Implemented, 
Equation  (2-7)  becomes: 

J 

i(t| )  =  X!  'x.(t|)  p.(t  )  (2-12) 

j  =  l  ^  ^  ^ 

Similarly,  Equation  (2-6)  describing  the  Pj^(t^)'s  becomes 

^z(ti  )|  i,Z(ti_i)  (Z-l'^j/Zl-l)  *  Pj(tl-l) 

Pj(ti)  =  -■  '  — - -  (2-13) 

J 

^z(ti  )  1  a,  Z(ti_i)  ^?-ll  Sk'^1-1  ^  P)c(^l-1^ 

and  Equation  (2-8)  similarly  is: 

( ti  )  I  a,  Z  ( ti  _i )  ^  ^1 1  i-j '  2.1-1 ) 

=  1/.;  (2n)n'/2  (Aj  (ti  )|l/2;  expC-(l/2)r5’(ti  )Ajl(ti  )rj  (ti  )3 

(2-14) 

and 

Aj(ti)  =  Hj(ti)P j(tX)H5(ti)  +  Rj(ti) 
rj(ti)  =  Zi  -  Hj(ti)£j(tX) 


m  =  the  dimension  of  z  (number  of  measurements) 

Rj  =  the  measurement  noise  covariance  in  the 
elemental  filter. 

2.3.2.  Sliding  the  Moving  Banlc  C3:253.  The 
moving  bank  MMAE  is  a  smaller  version  of  the  full-bank 
MMAE,  with  the  moving  bank  centered  around  a  parameter 
estimate.  Typically,  the  moving  bank  is  not  initially 
centered  on  the  true  parameter  point,  and  the  true 
parameter  point  may  change.  This  makes  the  decision  logic 
for  moving  the  "bank”  of  critical  Importance.  Several 
algorithms  have  previously  been  investigated  for  this 
purpose,  including  Residual  Monitoring,  Parameter  Position 
Estimate  Monitoring,  Parameter  Position  and  Velocity 
Estimate  Monitoring,  and  Probability  Monitoring  L3;  6;  19; 
9:  . 

2. 3, 2.1.  Residual  Monitoring.  Let  a  likelihood 
quotient  for  each  elemental  filter,  Lj(tj),  be  defined 
as  the  quadratic  form  appearing  in  Equation  (2-8): 

Lj(t^)  =  rj(t^  )A~^(t^  )r^ (t^  )  (2-15) 

The  decision  is  made  to  move  the  bank  if  at  time  tj^: 


minCL, ( t,  ) ,  L- ( t,  ) ,  .  .  . , 


L.(t, )D  >  T 


(2-16) 


where  T  Is  a  threshold  level  with  a  numerical  value  that  is 
determined  during  performance  evaluations.  If  the  true 
parameter  vector  value  is  outside  the  moving  bank,  it 
would  be  expected  that  all  the  likelihood  quotients  exceed 
the  threshold.  The  bank  is  then  moved  in  the  direction  of 
parameter  space  corresponding  to  the  filter  with  the 
smallest  L^,  as  that  filter  would  be  expected  to  be 
nearest  to  the  true  parameter  set.  This  method  should 
respond  quickly  to  a  real  need  to  move  the  bank  but  also 
gives  erroneous  results  for  a  single  instance  of  large 
residuals,  possibly  due  to  noise  corruption. 

2. 3. 2. 2.  Probability  Monitoring .  This  method  is 
similar  to  residual  monitoring  except  that  the  conditional 
hypothesis  probabilities,  generated  by  Equation  (2-6),  are 
monitored.  If  the  conditional  hypothesis  probability 
associated  with  an  elemental  filter  is  larger  than  a 
previously  determined  threshold,  the  bank  is  centered  on 
that  filter.  Maybeck  and  Hentz  found  this  decision  logic, 
as  well  as  parameter  position  monitoring,  to  provide  the 
best  performance.  However,  probability  monitoring  required 
fewer  computations  than  parameter  position  monitoring 
C6:93-99: 

2. 3. 2. 3.  Parameter  Position  Estimate  Monitoring . 


This  method  centers  the  bank  around  the  current  estimate 


of  the  true  parameter  set,  which  is  given  by: 


J 

'a'(t.  )  =  E  P^(t  )  (2-17) 

1  j  =  l  J  J  ^ 

where  J  is  the  number  of  filters  Implemented  in  the  moving 
bank.  Movement  is  initiated  when  the  bank  is  not  centered 
on  the  point  closest  to  the  current  true  parameter  set 
estimate  C  3 : 261  . 

2 . 3 . 2 . 4  .  Parameter  Position  and  Velocity 
Estimate  Monitoring.  This  method  estimates  the  velocity  of 
the  parameter  position  using  the  history  of  parameter 
position  estimates.  The  velocity  estimate  is  used  to 
estimate  the  position  of  the  parameter  set  at  the  next 
sample  time.  The  bank  is  centered  at  this  estimate  of  the 
future  parameter  point,  thereby  adding  "lead"  into  the 
positioning  of  the  bank  C241.  Maybeck  and  Hentz  found  this 
decision  logic  performed  worse  than  parameter  position 
estimate  monitoring  or  probability  monitoring  C6:85; 

19:231,  not  providing  much  desired  lead  but  causing  reduced 
stability  in  the  bank  location. 

2.3.3.  Bank  Contraction  and  Expansion .  The  filters 
in  the  moving  bank  model  do  not  necessarily  have  to  be  at 
adjacent  discretized  parameter  values;  see  Figure  2-1. 

This  may  decrease  the  accuracy  of  the  initial  estimate  but 
it  increases  the  probability  that  the  true  parameter  set 
lies  within  the  bank,  thereby  enhancing  Initial 
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convergence.  Maybeck  and  Hentz  found  that  parameter  acqui¬ 
sition  performance  can  be  Improved  by  starting  the  moving 
bank  with  a  coarse  discretization  so  that  the  entire 
parameter  value  range  lies  within  the  bank  and  then 
contracting  the  bank  Into  a  finer  discretization  when  the 
parameter  covariance  (Equation  (2-10))  drops  below  some 
selected  threshold  C3:28;  6:26;  19:251. 

Karnick  proposed  another  method  C91  that  may  improve 
acquisition.  This  method  monitors  the  probability  asso¬ 
ciated  with  a  "side"  of  the  bank;  see  Figure  2-2,  and  would 
be  calculated  as: 


Pside(ti)  =  -  (2-18) 

where  is  the  probability  distribution  function 

of  the  Kalman  filter  which  is  determined  using 

Equation  (2-14).  Several  possibilities  exist  for  threshold 
logic.  If  the  probability  associated  with  a  side  falls 
below  a  certain  threshold,  it  can  be  moved  toward  the 
center  of  the  movable  bank.  Conversely,  if  the  proba¬ 
bility  associated  with  a  side  rises  above  some  threshold, 
the  remaining  three  sides  are  moved  toward  the  center  of 
the  movable  bank.  A  third  possibility  is  moving  in  all 
four  sides  If  the  summed  probability  of  all  the  "side” 
filters  Is  below  some  threshold. 
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The  bank  may  also  need  to  be  expanded  if  the  true 


parameter  value  undergoes  a  jump  change  to  a  point  outside 
the  range  covered  by  the  bank.  This  jump  change  can  be 
detected  by  residual  or  probability  monitoring.  For 
residual  monitoring,  the  likelihood  ratios  for  all  the 
Implemented  filters  are  expected  to  be  large  and  to  exceed 
some  threshold.  For  probability  monitoring,  it  is  expected 
that  the  conditional  hypothesis  probabilities  be  "close"  in 
magnitude;  however,  performance  analysis  must  be  performed 
to  confirm  this.  The  subsequent  bank  contraction  is 
accomplished  in  the  same  manner  as  discussed  in  the 
previous  paragraphs. 

2.3.4.  Initialization  of  New  Elemental  Filters 
C3;29-31;  6:26-30;  9:29-321.  When  the  decision  is  made  to 
move,  expand,  or  contract  the  bank,  new  filters  must  be 
brought  on  line  and  "incorrect  filters"  discarded.  New 
filters  require  new  values  for  <t> ,  B^,  K  (steady  state 
Kalman  gain  matrix),  H,  ^(t|),  and  p^Ctj^).  Except  for  the 
last  two  terms,  these  are  predetermined  values  associated 
with  the  particular  filter  being  Implemented  because  for 
this  research  the  system  model  is  assumed  to  be  linear  and 
time  invariant. 

The  current  moving  bank  estimate  of  "^(t^)  is  an 
appropriate  choice  for  ^(t|)  for  a  new  elemental  filter. 
The  value  for  Pj(t^)  is  dependent  on  the  number  of  new 
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as  shown 


filters  being  implemented.  If  the  bank  "slides," 
in  Figure  l-3a,  this  involves  either  three  or  five  new 
filters.  The  probability  weighting  of  the  discarded 
filters  can  be  redistributed  equally  among  the  new  filters 
or  in  a  manner  that  indicates  the  estimated  "correctness" 
of  the  new  filter  and  still  maintains  the  sum  of  the 
probabilities  equal  to  one.  Hentz  suggested  the  following 
C7  :  293  : 

fj(z(ti))(l  - 

Pjch(ti)  =  - - -  (2-19) 

Z  fk(z(ti)) 

where  ch  =  changed,  unch  =  unchanged,  and  where  fj(s.(tj)) 
is  defined  in  Equation  (2-14)  but  with  the  residual 
replaced  by; 

rj(ti)  =  z^-  H^x^(tj)  (2-20) 

However,  this  requires  additional  computations  and  has 
demonstrated  no  significant  performance  Improvement  over 
dividing  the  probability  weighting  equally  among  the 
changed  filters  C6:1043. 

A  bank  expansion  or  contraction  can  result  in  the 
resetting  of  all  the  filters  in  the  bank  as  shown  in 


Figure  l-3b.  Dividing  the  probability  weighting  equally 
among  the  new  filters  is  appropriate  since  the  old 
^  probability  weightings  may  no  longer  be  valid. 

i 

2.4.  Controller  and  Estimator  Design  13:18-22;  6:33-43; 

9  :  32-393 

Several  controller  and  estimator  designs  are 
appropriate  for  Implementation  with  the  moving  bank  or 
full-bank  MMAE .  All  designs  considered  use  the  "assumed 
certainty  equivalence  design"  technique  123:2411,  which 
consists  of  developing  an  estimator  cascaded  with  a  deter¬ 
ministic  full-state  feedback  optimal  controller.  This 
method  assumes  Independence  between  controller  and  esti¬ 
mator  design  and  is  the  optimal  stochastic  controller 
design  for  a  linear  system  driven  by  white  Gaussian  noise 
with  quadratic  performance  criterion  123:173. 

The  moving-bank  MMAE  or  the  fixed  bank  MMAE  is  the 
estimator  used  in  this  thesis.  Each  elemental  estimator 
within  the  full  bank  Is  a  constant-gain  Kalman  filter 
associated  with  a  particular  point  in  the  parameter  space. 
Each  design  assumes  a  time  invariant  system  with  stationary 
noise.  Propagation  of  the  elemental  filter  estimate, 

X  ( t ) ,  is  given  by: 

-^(tT)  =  «i>^(tr , )  +  B^^u(t,_, ) 


{  2-21  ) 


and  the  estimate  is  updated  by: 


x^(t|)  =  x^(t^)  +  Kj^[:z(t^)  -  H^x^(t*):  (2-22) 

where  the  subscript  "k"  Indicates  association  with  a 
particular  point  a^^  In  the  parameter  space. 

The  design  of  each  controller  Is  similar.  Each  is 
a  linear,  quadratic  cost,  (LQ)  full-state  feedback  optimal 
deterministic  controller,  based  on  a  time  invariant  error 
state  space  formulation  and  constant  cost  weighting  matri¬ 
ces  in  the  quadratic  cost  to  be  minimized  C21:2971.  The 
form  of  the  cost  equation  to  be  minimized  is  as  follows: 

^  T 

J  =  E-  X)  1/2CX  (t. )Xx{t, )  +  u  (t  )Uu(t. )1  + 
i  =  0 

l/2x'^(t^^j^  )XjX(t^^j^)[  (2-23) 

where  : 

J  =  cost  function  to  be  minimized 
X  =  state  weight  matrix 
X^  =  final  state  weight  matrix 
U  =  control  weight  matrix 

N  =  number  of  sample  periods  from  t.,  to  t^ 


» 


=  final  time 

=  last  time  a  control  is  applied  and  held  constant 
over  the  next  sample  period 


The  controller  is  steady-state  constant-gain,  with  gains 
dependent  upon  the  particular  value  of  the  parameter  set 
used  in  the  design.  The  LQ  controller  is  developed  fully 
in  Appendix  A. 

Three  estimator/controller  combinations  have  been 

researched  in  past  theses.  First,  the  estimator  provides 

only  a  state  vector  estimate  to  a  fixed-gain  controller 

which  is  designed  around  a  nominal  value  of  the  uncertain 

parameter  set,  a  .  The  controller  algorithm  is  of  the 

nom 

f  orm : 


u(t,  )  =  -G*Ca  :  ) 

1  c  ^om  1 


( 2-24  ) 


The  second  design  method  is  for  the  estimator  to 
provide  parameter  and  state  vector  estimates  to  a 
controller  with  gains  that  are  dependent  on  the  parameter 
estimate : 


u(t^)  = 


-G^Ca(t,  . )]  x(t, ) 
c  1-1  —  i 


(2-25) 


tSiil, 


I 


where  the  parameter  estimate  generated  at  the  previous 
sample  time  is  used  in  order  to  reduce  computational  delay. 

The  third  approach  is  to  implement  an  elemental 
controller  for  each  of  the  elemental  filters  of  the  sliding 
bank.  The  control  outputs  are  probabilistically  weighted, 
similar  to  Equation  (2-12),  to  form: 

J 

u(t.)  =  (2-26) 

1  =  1 

where  : 

Uj(ti)  =  -G*Cajl  x^(t|)  (2-27) 

This  is  usually  referred  to  as  a  multiple  model  adaptive 
controller  (MMAC)  [23:2533.  This  research  uses  the  last 
approach.  However,  initially  the  controller  design  shown 
in  Equation  (2-27)  was  also  going  to  be  Implemented  in  this 
research.  The  nominal  controller  was  needed  because  If 
adaptive  control  Inputs  were  supplied  to  the  system  before 
an  accurate  estimate  of  the  true  parameter  was  obtained, 
the  control  inputs  could  drive  the  system  unstable. 
Therefore,  a  nonadaptive  controller  was  designed  that  was 
very  robust  to  parameter  variation  at  the  expense  of 
controller  performance  so  that  the  controller  Inputs 


would  at  least  maintain  stability  regardless  of  the 
estimate  of  the  true  parameter. 


2.5.  Ambiguity  Funct Ion  Analysis  r3;  6;  9;  20:97-991 

Ambiguity  function  analysis  can  provide  information 
about  the  performance  of  an  estimator.  The  generalized 
ambiguity  function  is  given  by: 


A^ (a,at  ) 


LCa^Z^I  f z ( t i )  I  a ( t i ) ( ?i I ?t ^l^i 

( 2-28  ) 


where  a  is  the  parameter  vector,  a^  is  the  true  parameter 

vector,  and  LCa,Z^l  Is  a  likelihood  function  upon  which  a 

parameter  estimate  would  be  based  via  maximum  likelihood 
techniques.  For  a  given  value  of  a^,  the  curvature  of  the 

function  of  a,  at  the  value  of  a^,  provides  information  on 

the  ability  of  the  filter  to  estimate  that  parameter:  the 

sharper  the  curvature,  the  greater  the  precision.  This 

curvature  is  Inversely  related  to  the  Cramer-Rao  lower 

bound  on  the  estimate  error  covariance  matrix  by 


Et(a  -  at)(a  -  at)"^!  >  C-{  d^/^a^)  A^  ( a ,  at  )  I  a=a^  ^  ‘  ^ 

■  '  12-29) 

The  ambiguity  function  value  A^(a,a^)  for  any  a  and 
a^  can  be  calculated  from  the  output  of  a  conventional 
nonadaptlve  Kalman  filter  sensitivity  analysis  122:97-991 
in  which  the  "truth  model"  Is  Identical  to  the  model  upon 


which  the  Kalman  filter  is  based,  except  that  they  are 
based  on  and  a,  respectively.  The  ambiguity  function  is 
then  given  by 

A”(a,^)  =  52  {{-m/2  ln(2n)  -  1/2  IntlACt  ;a)|: 

^  j  =  i-M+l\ 

-  1/2  tr  |A"^(tj;a)[H(tj  )Pg(t^;a^,a)H'^(tj  )  +  R(tj): 

-  n/2  ln(2n)  -  1/2  InC I P ( t[ ; a  )  I  1 

-  1/2  trCP'^  (t!';a  )Pg  (t|;a^,a  )3  (2-30) 

where  : 

A(ti,a)  =  EH(ti  )P(ti;;a)H'r(ti  )  +  R(ti)l"l 

for  the  Kalman  filter  based  on  a, 

Pg(t|^;at,a)  =  the  covariance  matrix  of  the  error 

between  the  state  estimate  of  the  Kalman 
filter  based  on  a  and  the  states  of  the 
true  system  based  on  a^,  where  "  "  or  "+" 
denotes  before  or  after  Incorporation  of  the 
1^^  measurement 

m  =  the  number  of  measurements 

n  =  the  number  of  states 

M  =  the  number  of  the  most  recent  residuals 

included  in  the  likelihood  function  (rather 
than  using  all  residuals  back  to  the  first 
sample  term 

The  terms  are  summed  over  the  most  recent  M  sample  times 
[22:983;  however,  here  M  is  set  equal  to  one.  This 
approximation  can  be  made  because  the  summation  term  is  the 
component  of  Equation  (2-30)  that  is  sensitive  to 


parameter  variation,  while  the  rest  of  Equation  (2-30)  is 
sensitive  to  the  state  estimation.  Therefore,  the  general 
shape  of  the  curve  can  be  found  by  using  only  the  first 
term  in  the  summation  term  (see  Equation  2-31).  This 
reduces  the  size  of  the  fluctuations  in  the  value  of 
Consequently,  this  flattens  the  surface  of 
the  plot  of  the  ambiguity  function  plotted  as  a  function 
of  the  parameter  a.  The  main  benefit  of  setting  M  =  1  is 
that  this  significantly  reduces  the  number  of  computations. 

Flllos  encountered  numerical  difficulties  while 
evaluating  the  ambiguity  function.  The  numerical 
difficulties  were  caused  by  numerical  precision  problems 
encounter  when  cascaded  multiplies  were  used  In  determinant 
evaluation  i:3:64l.  Therefore,  it  Is  Impossible  to  compute 
the  ambiguity  function  as  described  In  Equation  (2-30). 

The  numerical  difficulties  are  overcome  by  approximating 
the  expressions  for  the  probability  weighting  factors 
(Equation  2-14)  and  the  ambiguity  function  (Equation  2-30). 
Equation  (2-14)  is  approximated  as 

fj(z(ti))  =  expC-(l/2)r5’(ti  )Ajl(ti  )r  (ti  )3 

(2-31) 

This  is  no  longer  a  true  density  function  because  the 


scale  factor  is  Incorrect:  however,  because  of  the  deno¬ 
minator  in  Equation  (2-13),  the  probability  weightings  are 
still  correct  in  the  sense  that  they  add  to  one  C3:65].  If 
the  determinants  of  the  A  matrices  of  the  elemental  filters 
are  expected  to  be  approximately  equal  in  magnitude,  in  the 
absence  of  numerical  problems,  the  relative  magnitudes  of 
the  value  of  the  ambiguity  functions  are  not  significantly 
altered  13:651.  Equation  (2-30)  is  approximated  by 
removing  the  terms  containing  the  determinants  of  P(t^) 
and  A(t^).  Equation  (2-30)  becomes 

Aj(a,a^.)  =  -m/2  ln(2n)  -  n/2  ln(2n) 

-1/2  trCA"^(t.  ;a)CH(t.  )Pg  (t7;a^,a)H'^(t.  )  +  R(t.)33 
-1/2  trCP"^(t^;a)Pg(t[;a^,a)3  (2-32) 

However,  if  M  is  equal  to  one,  then  the  last  term  is 
ignored  and  Equation  (2-32)  becomes: 

Aj{a,aj.)  =  -m/2  ln{2n) 

-1/2  trCA"^(t^;a)CH(t^  )Pg(t';a^,a)H'^(t^  )  +  R(t^)]3 

(2-33) 

This  is  a  reasonable  approximation  since  the  determinants 
of  P(t|)  and  A(t^)  have  minimal  effect  on  the  ambiguity 
function.  Therefore,  its  primary  sensitivity  is  in  the 


functions  that  are  being  preserved  r3:663.  Karnick  used 
this  approximation  in  his  thesis  and  the  ambiguity  plots  he 
produced  contained  numerous  peaks  and  valleys,  making  it 
very  difficult  to  determine  the  correct  space 
d iscret i zat i on . 

However,  there  are  other  techniques  that  can  be  used 
to  determine  the  discretization  level.  First,  Insight  can 
be  gained  from  the  values  of  the  dynamics  noise  strength 
(Q)  and  the  measurement  noise  covariance  (R).  Large  Q  and 
R  values  indicate  a  need  for  large  space  discretization 
because  such  large  Q  and  R  values  mask  differences  in  the 
filters  in  a  tightly  spaced  discretization.  Another 
technique  is  to  conduct  a  sensitivity  study  by  giving  the 
filter  an  incorrect  parameter  value  and  then  observing  the 
state  estimate.  The  more  Inaccurate  the  state  estimate 
becomes,  for  a  given  size  variation  between  filter-assumed 
and  true  parameter  values,  the  tighter  the  space 
discretization  needs  to  be  for  that  parameter. 


2.6.  Summary 

This  chapter  developes  algorithms  necessary  for 
implementation  of  the  full-scale  and  moving-bank  multiple 
model  adaptive  estimator  and  an  appropriate  adaptive 
controller  based  on  this  type  of  estimation.  The  moving 
bank  MMAE  is  expected  to  yield  significant  computational 
savings  over  the  full-scale  MMAE,  if  they  are  both  based 


on  the  same  discretization  level.  However,  a  viable 
alternative  to  a  moving  bank,  which  is  also  investigated, 
is  a  fixed-bank  algorithm  based  upon  a  coarser  discreti¬ 
zation.  The  ambiguity  function  analysis  is  also  developed. 

Ambiguity  functions  are  expected  to  give  insight  into  the 
parameters  that  need  adaptive  estimation  and  into  the 
appropriate  levels  of  discretization  of  the  parameter  space 
to  perform  such  estimation. 


III.  Rotat inq  Two-Bav  Truss  Model 

3.1.  Introduction 

This  chapter  develops  the  system  equations  for  a 
large  flexible  space  structure  modeled  as  a  two-bay  truss. 
The  truss  is  allowed  to  rotate  around  a  fixed  point, 
thereby  incorporating  both  rigid  body  rotation  and  bending 
mode  dynamics.  The  system  model  is  first  developed  in 
physical  variables  and  then  transformed  into  the  modal 
variable  format.  This  is  followed  by  a  complete  physical 
description  of  the  two-bay  truss.  This  model  is  taken  from 
Karnick's  thesis  C93.  Section  five  motivates  the  need  for 
order  reduction  and  develops  the  order  reduction  technique 
used  in  this  thesis. 


3.2.  Second  Order  and  State  Space  Form  Models 

The  forced  vibration  of  a  large  space  structure  with 
active  controls  and  n  frequency  modes  can  be  described  with 
the  following  general  second-order  differential 
equations  C17;  321: 


Mr(t)  +  Cr(t)  +  Kr(t)  =  Fi(u,t)  +  F2(t)  (3-1) 

where  : 

M  =  constant  nxn  mass  matrix 
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C  =  constant  nxn  damping  matrix 


K  =  constant  nxn  stiffness  matrix 

r.(t)  =  n-vector  representing  structure's  physical 
coord i nates 

Fi(u,t)  =  control  input 
F 2 ( t )  =  disturbances 


The  control  system  is  assumed  to  consist  of  a  set  of 
discrete  actuators.  The  location  of  the  actuators  are 
given  in  Section  3.4.4.  The  external  disturbances  are 
represented  by  white  noise,  thus  producing: 


Mr ( t )  +  C r ( t )  +  Kr(t)  =  -bu(t)  -  gw  (3-2) 


where  denotes  a  vector  stochastic  process  and: 

u ( t )  =  vector  of  dimension  m  representing  actuator 
i  nput 

b  =  nxm  matrix  identifying  position  and  velocity 
relationship  between  actuators  and  controlled 
var iables  r 17 1 

w  =  vector  of  dimension  r  representing  dynamics 

~  driving  noise,  where  r  is  the  number  of  noise 
inputs 

g  =  nxr  matrix  identifying  position  and  rela¬ 
tionships  between  dynamics  driving  noise  and 
controlled  variables. 


The  state  representation  of  Equation  (3-2)  can  be 
written  as: 
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X  =  FX  +  Bu  +  Gw 


(3-3) 


where  the  2n-dimensional  vectors  x  and  x  are  given  by: 


(  3-4  ) 


and  the  open-loop  plant  matrix  F,  the  control  matrix  B,  and 
the  noise  matrix  G  are  given  by: 


-M  -M 


B  =  0 


2nx2n 


2nxm 


(  3-5) 


G  =  r  0 


-M  ^g 


2nxr 


It  is  assumed  that  the  noise  enters  the  system  at  the  same 
place  as  the  actuators;  therefore,  the  b  matrix  is  equal  to 
the  g  matrix  and  m  is  equal  to  r.  In  addition,  it  is 
assumed  measurements  are  available  from  accelerometers  and 
gyros.  The  location  of  the  sensors  are  given  in  Section 
3.4.4.  The  output  signals  of  the  accelerometers  are 
integrated  once  to  provide  velocity  measurements  and  once 
again  to  provide  position  measurements.  Therefore,  the 
position  and  velocity  measurements  are  assumed  co-located. 


These  position  and  velocity  measurements  are  then  fed  into 
the  system  model.  It  is  also  assumed  that  the  roundoff 
errors  due  to  the  A/D  conversion  associated  with  the  two 
integrations  are  small  and  can  be  accounted  for  by 
increasing  the  measurement  noise  (R).  In  addition,  it  is 
assumed  that  the  position  and  velocity  measurements  are 
noise  corrupted  due  either  to  deficiencies  in  the  model  of 
the  sensor  or  some  actual  external  measurement  noise.  The 
position  and  velocity  measurements  are  modeled  as: 

H  0  X  +  y  { 3-6 ) 

•*««/ 

0  H ' 

px2n 

where : 

p  =  the  number  of  measurements 
V  =  an  uncertain  measurement  disturbance  of 

dimension  p  and  modeled  as  a  white  Gaussian 
noise  C21:114] 

H  =  position  measurement  matrix 

* 

H  =  velocity  measurement  matrix 

Since  the  position  and  velocity  measurements  are 

t 

co-located,  the  measurement  matrices  are  equal  (H  =  H  ). 


3.3. 


Modal  Analysis 

Modal  analysis  is  used  to  transform  the  system  into  a 
set  of  independent  equations  by  transforming  the  system 
from  physical  coordinates  to  modal  coordinates.  In  order 
to  achieve  decoupling,  the  damping  matrix  must  be  assumed 
to  be  a  linear  combination  of  the  mass  and  stiffness 
matrices  C16]: 

C  =  a  K  +  ^ M  ( 3-7 ) 

The  determination  of  a  and  $  is  discussed  later  in 
this  section. 

The  modal  coordinates  are  related  to  the  physical 
coordinates  by 

r  =  T  n  (  3-8  ) 

where  r.  is  as  defined  previously  and  Q  represents  the 
modal  coordinates.  The  matrix,  T,  is  an  nxn  matrix  of 
eigenvectors  and  is  the  solution  to  C17;  27;  32;  333: 

2 

w  MT  =  KT  ( 3-9  ) 

The  values  of  w  which  satisfy  Equation  (3-9)  are  natural  or 
modal  frequencies  18:663.  Substituting  Equation  (3-8)  into 
Equation  (3-3)  gives 


•  I  It  t  I 

X  =  Fx  +  Bu  +  Gw 


( 3-10) 


where  x  Is  now  defined  as: 


2nxl 


Q 

h' 


( 3-11  ) 


2nxl 


and  the  open  loop  plant  matrix  F  ,  the  control  matrix  B  , 

I 

and  the  noise  matrix  G  are: 


r 

0 

I 

■ 

-1 , 

-1  -1 

-T  M  KT 

-T  M 

CT 

2nx2n 

•  «« 

0 

1 

G  = 

- 

0 

^-1  -1^ 

-1  -1 

L -T  M  b  _ 

T  M  b 

(3-12) 


2nxm 


2nxm 


The  F  matrix  can  also  be  written  in  the  following 
form  C32;  3  33  : 


I 

F 


0 


I 


3 


C-2fi«j3 


(3-13) 


2nx2n 


where  each  of  the  four  partitions  are  nxn  dimensional  and 
diagonal  and  where  and  are  the  undamped  natural 
frequency  and  the  damping  ratio  of  the  ith  mode.  The  mea¬ 
surements  become: 

I 

X  +  V  ( 3-17 ) 

Since  the  system  is  formulated  in  modal  coordinates 
C171,  it  is  assumed  that  uniform  damping  exists  throughout 
the  structure.  The  level  of  structural  damping  is  deter¬ 
mined  by  selecting  a  value  for  the  damping  coefficients 
(  ^ )  and  substituting  this  value  into  Equation  (3-13). 

The  particular  value  of  the  damping  coefficient  has  no 
effect  on  the  calculation  of  since  it  is  the  natural 
or  modal  frequency.  The  assumption  simplifies  the 
determination  of  structural  damping  and  allows  a  better 
physical  insight  into  formulating  the  problem  than  does  the 
selection  of  values  for  a  and  0  as  shown  in  Equation 
(3-7).  The  damping  coefficient  of  =  0.005  is  chosen 
for  Implementation  because  it  is  characteristic  of  damp¬ 
ing  associated  with  large  space  structures  C17;  241. 

3.4.  Two -Bay  Truss 

3.4.1.  Introduction .  This  section  describes  the 


physical  structure  of  the  two-bay  truss  rotating  about  a 


fixed  point.  This  section  discusses  the  physical  dimen¬ 
sions  of  the  model,  the  analysis  used  to  develop  the 
two-bay  truss  model,  the  sensors  and  actuators,  and  the 
physical  parameter  variations  to  be  considered  in  the 
mode  1  . 

3.4.2  Background .  For  this  research  a  fixed  two-bay 
truss,  which  was  used  to  research  active  control  laws  for 
vibration  damping  fl73  (see  Figure  3-1),  was  modified  twice 
to  produce  the  system  model.  First,  non-structural  masses 
were  added  to  nodes  1  through  4,  which  reduced  the  struc¬ 
tural  frequencies,  thereby  making  the  problem  more  like  a 
large  space  structure  C171.  Second,  rigid  body  motion  was 
added  to  the  truss  (see  Figure  3-2).  These  modifications 
produced  a  model  which  approximated  a  space  structure  that 
has  a  hub  with  appendages  extending  from  the  structure. 

The  mass  of  the  hub  was  large  relative  to  the  mass  of  the 
appendage  and  the  hub  can  then  be  rotated  to  point  the 
appendage  in  a  commanded  direction. 

The  rigid  body  motion  was  established  by  adding  a 
point  (node  7)  that  remained  fixed  while  the  two-bay  truss 
was  free  to  rotate  about  this  point  in  the  x-y  plane;  see 
Figure  3-2.  The  truss  was  connected  to  this  point  with 
rods  having  radii  that  were  large  relative  to  the  rods  used 
to  construct  the  two-bay  truss.  Making  the  rods  large  made 
a  very  "stiff"  link  between  the  truss  and  node  7,  which 


introduced  high  frequency  modes  into  the  structure  but  kept 
the  lower  modal  frequencies  similar  to  the  case  where  the 
truss  was  fixed. 


3.4.3.  Two-Bay  Truss  Construct i on .  The  structure 
consisted  of  13  rods  which  were  assumed  to  be  constructed 
of  aluminum,  having  a  modulus  of  elasticity  of  10  psi  and 
weight  density  of  0.1  Ib/in  C323.  The  cross-sectional  areas 
of  each  member  shown  in  Figure  3-2  are  given  in  Table  3-1. 


TABLE  3-1 


STRUCTURAL  MEMBER'S 

CROSS- 

SECTIONAL  AREAS 

Member 

Area  ( i n^ ) 

Member 

Area  (in^) 

a 

.  00321 

h 

.00328 

b 

.00100 

i 

.00439 

c 

.00321 

j 

.00439 

d 

.01049 

k 

. 20000 

e 

.  00100 

1 

. 20000 

f 

.01049 

m 

. 20000 

g 

.  00328 

- 

- 

The  cross-sectional  areas  of  rods  a  through  j  were  cal¬ 
culated  by  optimizing  the  weight  of  the  structure  shown  in 
Figure  3-1.  The  large  space  structure  was  calculated  with 
its  weight  minimized  subject  to  the  constraint  that  the 
fundamental  frequency  remain  unchanged.  This  constraint 
was  needed  so  that  the  natural  frequency  response  of  the 
system  remained  constant.  Rods  k  through  m  were  used  to 


make  the  "stiff"  link  between  node  7  and  the  two-bay  truss. 
In  order  to  achieve  this  stiffness,  the  areas  of  rods  k  and 
1  were  arbitrarily  selected  to  be  large  relative  to  the 
area  of  other  rods. 

2 

Non-str uctura 1  masses  of  1.294  lb-sec  /in,  were 
located  at  positions  1,  2,  3,  and  4  as  shown  in  Figure  3-2. 
The  non-structural  mass  was  very  large  compared  to  the 
structural  mass  to  achieve  the  low  frequencies  associated 
with  large  space  structures  tl73.  The  actual  value  of  the 
non-structural  mass  was  selected  using  an  optimization 
technique  1333,  which  found  the  mass  necessary  to  attain  a 
frequency  of  0.5  Hz  in  the  lowest  mode  for  the  fixed 
two-bay  truss  C173. 

The  mass  and  stiffness  matrices,  describing  the 
system  model,  were  obtained  using  finite  element  analysis 
[313,  which  models  a  structure  as  consisting  of  a  finite 
number  of  nodes  connected  by  elements.  The  software 
program  that  was  used  C333  has  the  capability  to  use  a 
number  of  different  elements.  In  this  research  effort, 
the  rods  were  described  by  cross-sectional  area,  modulus  of 
elasticity,  and  weight  density.  The  finite  element  program 
produced  mass  and  stiffness  matrices  with  dimension  equal 
to  the  number  of  degrees  of  freedom  (DOF)  associated  with 
the  model.  Each  row  of  the  mass  and  stiffness  matrix  was 
associated  with  a  specific  node  and  DOF.  For  the  two-bay 


truss  shown  in  Figure  3-2,  row  1  of  each  mass  and  stiffness 
matrix  is  associated  with  the  x-axis  DOF  of  node  1.  Each 
node  actually  has  three  translational  DOF.  However,  only 
planar  motion  was  being  considered  (in  order  to  simplify 
the  model);  therefore,  the  nodes  were  modeled  with  only  two 
DOF.  For  this  problem,  node  7  was  fixed.  Therefore,  all 
three  DOF  associated  with  this  node  were  eliminated, 
thereby  reducing  the  dimensionality  of  the  mass  and 
stiffness  matrices  to  12  states  yielding  a  24-state  model. 
The  mass  and  stiffness  matrices  for  the  specifications 
previously  discussed,  are  listed  in  Appendix  A.  These  were 
the  nominal  matrices  from  which  parameter  variations  were 
cons  idered  . 

3.4.4.  Sensors  and  Actuators .  Accelerometers 
were  located  at  nodes  1  and  2  as  shown  in  Figure  3-2.  The 
accelerometers  should  not  be  located  at  a  node  of  the 
bending  modes  being  detected  because  the  displacement  of 
the  truss  caused  by  the  bending  modes  can  not  be  detected. 
The  accelerometer  outputs  were  integrated  twice  to  produce 
velocity  and  position  measurements.  There  were  also  two 
additional  gyro  sensors  co-located  on  the  hub  (node  7)  of 
the  two-bay  truss.  These  sensors  provided  angular 
displacement  and  velocity.  Actuators  were  placed  at  nodes 
1  and  2  as  shown  in  Figure  3-2.  The  actuators  do  not  have 
to  be  co-located  with  the  sensors,  however  they  are  in  this 


case  to  simplify  the  system  model.  These  actuators  were 
assumed  to  be  force  sensors  (thrusters,  etc.).  An 
additional  actuator  was  located  on  the  hub  and  is  assumed 
to  be  an  inertia  wheel. 

The  states  corresponding  to  velocity  and  position 
were  directly  available  in  physical  variable  formulation, 
(Equations  3-3,  3-4,  3-5)  while  the  states  corresponding 
to  angular  displacement  and  velocity  were  directly 
available  in  modal  formulation  (Equations  3-10,  3-11,  3-12, 
3-14).  The  H  and  b  matrices  were  constructed  by  cal¬ 
culating  separate  matrices  in  the  different  state  space 
formulations.  These  matrices  were  augmented  after  the 
physical  variable  formulations  have  been  transformed  into 
modal  coordinates. 

3.4.5.  Physical  System  Parameter  Uncertainty.  The 
purpose  of  this  thesis  was  to  test  the  moving-bank  multiple 
model  adaptive  estimation  and  control  algorithms.  There¬ 
fore,  the  model  must  have  parameter  uncertainty  which 
allowed  adaptive  estimation  to  be  applied.  A  10  by  10 
point  parameter  space  was  created  by  considering  two 
physically  motivated  parameter  variations.  The  parameters 
that  were  varied  were  the  masses  of  the  four  non-structural 
masses  and  the  stiffness  matrix.  The  weight  variation  can 
be  physically  related  to  fuel  being  expended  from  or  added 
to  a  tank  or  weight  being  shifted  to  a  different  section 


(other  than  the  two-bay  truss)  of  the  space  structure.  The 
stiffness  variation  can  be  associated  with  structural 
fatique  in  the  rods  or  a  failure  of  a  member  within  the 
structure  itself.  The  exact  parameter  discretization  will 
be  determined  in  a  sensitivity  study  performed  in  this 
thesis.  Initially  in  forming  the  parameter  space,  the  four 
non-structural  mass  were  multiplied  by  constant  terms  that 
ranged  from  0.5  to  1.4,  while  the  stiffness  matrix  was 
multiplied  by  constant  terms  that  ranged  from  0.8  to  1.16. 
However  following  the  space  discretization  study,  the 
masses  were  multiplied  by  constant  terms  that  ranged  from 
0.5  to  1.5,  while  the  stiffness  matrix  was  multiplied  by 
constant  terms  that  ranged  from  0.5  to  1.4.  The  mass  and 
stiffness  matrices  were  chosen  as  parameters  because  they 
appeared  in  the  system  equations  and  made  logical  sense; 
however,  with  further  study  more  appropriate  parameters  may 
be  determined.  An  additional  technique  used  in  determining 
the  appropriate  uncertain  parameters  is  ambiguity  function 
plots  n  33 . 

3.5.  State  Reduct  1  on 

3.5.1.  Introduction.  The  mass  and  stiffness 
matrices  were  previously  shown  to  be  of  dimension  12.  This 
produces  a  system  model  that  has  24  states,  which  was  much 
larger  than  desired  for  the  control  synthesis  of  this 
thesis  effort  and  for  a  practical  control  application. 


This  section  develops  a  method  of  order  reduction  referred 
to  as  singular  perturbations  CIO;  11;  17;  23:2191,  which 
assumes  that  faster  modes  reach  steady  state  essentially 
instantaneously.  This  section  discrlbes  the  method  of 
singular  perturbations  and  then  discusses  the  magnitude  of 
the  order  reduction. 

3.5.2  Deve 1 opment .  The  deterministic  system  is 
reformulated  as  follows: 


(3-15) 


(3-16  ) 


The  states  are  to  be  retained  and  and  F22 

are  square  matrices.  If  the  high  frequency  modes, 
represented  by  states,  are  eliminated,  steady  state  is 
assumed  to  be  reached  instantaneously  for  these  modes 
( X 2  =  0.)-  The  X2  states  are  then  expressed  in  terms 
of  the  x^  states  and  controls,  where  it  is  assumed  F  ^ 
exists,  as  follows: 


?2  “  -  ~  ^21-1^  ^22-2  ^  ®2-  (3-17) 

-  =  -F.,^(F-.,x  +  B_u) 


X 


(3-18) 


Substituting  for  X2  gives 


^^1 

II 

IX 

+  Bu  z  =  Hx,  +  Du 

r-  -  r-1  r- 

(3-19) 

where ; 

F 

r 

II 

-  F  F~^F  ) 

122221 ' 

(3-20) 

B 

r 

=  (B^ 

^12^22^2* 

(3-21) 

H 

r 

=  (H, 

”2^22^21’ 

(3-22) 

D 

r 

=  (-H^ 

^22®2  ’ 

(3-23) 

Note  that  the  matrix  did  not  exist  before  order 
reduction.  It  is  a  direct-feed  forward  term  which  is  not 
in  the  unreduced  system  C173. 

This  order  reduction  technique  is  now  applied  to  a 
system  of  the  form  of  Equation  (3-10).  Reordering  Equation 
(3-13)  into  the  reduced-order  form  produces  Equation 
(3-24),  where  the  upper  partition  contains  the  modes  to  be 
retained  while  the  lower  partition  contains  those  assumed 
to  reach  steady  state  instantaneously: 


(  3-24 ) 

2nx2n 

The  exact  dimension  of  is  equal  to  two  times  the  number 
of  modes  contained  in  the  system  model.  The  exact 
dimension  of  is  equal  to  two  times  the  number  of  modes 
assumed  to  reach  steady  state.  Comparing  Equation  (3-24) 
to  Equation  (3-15)  shows  that  the  partitions  ^21 

are  zero.  Substituting  this  result  into  Equation  (3-19) 
yields  ; 


F  = 


t-2f2-2^ 


(  3-25) 
(  3-26  ) 


Hj.  =  (3-27) 

Dr  =  (-H2F22B2)  (3-28) 


is  the  only  term  in  Equation  (3-19)  to  Equation  (3-24) 
that  is  dependent  upon  terms  associated  with  the  states 
assumed  to  reach  steady  state  instantaneously.  The  other 
reduced-order  matrices  are  calculated  simply  by  truncating 


those  states  associated  with  x-. 


Calculation  of  can  be  greatly  simplified  by 
examining  the  form  of  Equation  (3-28)  H2  is  similar  in 
form  to  Equation  (3-6): 


(  3-29  ) 


represents  measurement  of  the  unmodeled  position  states 

t 

while  represents  measurement  of  the  unmodeled  velocity 
states.  In  Equation  (3-6),  it  is  assumed  that  the  posi¬ 
tion  and  velocity  measurement  matrices  are  identical 
because  the  measurments  are  derived  from  an  acceleration 
measurement.  The  same  assumption  can  be  made  in  Equation 
(3-29);  however,  the  distinction  between  the  velocity  and 
measurement  matrices  is  retained  since  it  is  shown  in 
Equation  (3-33)  to  be  important  in  the  general  development 
of  the  reduced  order  matrices.  As  is  shown  in  Equation 
(3-27),  F22  is  3  square  matrix  of  the  form: 


F 

22 


0 


C-2 


(3-30) 


where  each  of  the  four  partitions  is  a  square,  diagonal 
matrix  whose  dimension  is  dependent  upon  the  number  of 


states  to  be  retained.  Its  inverse  is  C83: 


^ ,  «•  .  «•'_  W-  C-  ■"-  ^ 


,-l 


22 


C-u;2]-1[2  <-2^2’ 


l-w^J  ^ 


I 


(  3-3]  ) 


B^  is  similar  in  form  to  the  matrix  B  described  in  Equation 
(  3-5  )  : 


0 

b 


(3-32) 


where  b  represents  the  rows  of  the  matrix  product  -M  "^b 
corresponding  to  the  unmodeled  states.  Evaluation  of 
Equation  (3-28)  yields: 


O'  2,-1^' 

He  -2  b 


0 


( 3-33  ) 


pxm 


where  p  is  the  number  of  measurements.  Only  the  position 
measurements  are  affected  since  the  lower  portion  is  zero. 
The  matrix  is  only  dependent  upon  the  position  portion 
of  the  measurement  matrix  and  not  the  velocity  measurement 


matrix.  The  inverse  of  is  easily  calculated  since 


the  matrix  is  diagonal.  An  example  of  detailed  system 
matrix  development  and  order  reduction  is  listed  in 


Appendix  B. 


3.5.3.  Order  Reduction  Selection .  The  number  of 
modes  retained  was  determined  by  examination  of  the 
eigenvalues  and  frequencies  of  the  unreduced  system 
(Table  3-2).  The  eigenvalues  were  determined  for  an 
undamped  system.  This  greatly  simplified  the  problem  and 
was  a  good  approximation  because  large  space  structures 
have  very  small  damping  coefficients  (  (^=  0.005).  The 
frequencies  can  be  distinctly  divided  into  several  groups 
of  closely  spaced  frequencies.  For  example,  modes  4,  5, 
and  6  were  clearly  one  set  of  closely  spaced  frequencies. 
When  reducing  the  order  of  system  by  the  method  of  singular 
perturbations,  it  is  not  desired  to  make  the  reduction  at  a 
point  which  divides  a  group  of  "closely  spaced"  frequen¬ 
cies  r241.  At  the  same  time,  a  sufficient  number  of  fre¬ 
quencies  must  be  retained  in  order  to  do  an  adequate  job  of 
estimation  and  control.  In  this  research,  the  truth  models 
for  simulating  the  real  world  and  design  models  upon  which 
to  base  the  est 1 mator /control  let  algorithms  were  inves¬ 
tigated  at  dimensions  of  6  for  the  truth  models  and  6  for 
the  filter  models.  Therefore,  this  research  did  not 
analyze  the  effect  of  order  mismatch  between  the  filter 
models  and  the  truth  model. 

3.6.  Summary 

This  chapter  developed  the  system  equations  for  the 


two-bay  truss  with  rigid  body  motion.  The  mathematical 


Mode  No. 


Eigenvalues 


Frequencies  (Hz) 


1 

0,0000 

0.0000 

2 

8.8922 

1.4152 

3 

22.5492 

3.5888 

4 

29 . 5444 

4.7021 

5 

31.1519 

4.9580 

6 

32.8002 

5.2203 

7 

54.3893 

8.6563 

e 

58.1592 

9.2563 

9 

985.9204 

156.9141 

10 

9018.8987 

1435.4023 

11 

11515.9941 

1832.8274 

12 

19956 . 5072 

3176.1768 

model  was  dependent  upon  physical  parameters  which,  In 
reality,  vary  from  those  used  in  the  mathematical  model. 
The  moving-bank  MMAE  was  used  to  estimate  both  the  reduced 
order  system  states  and  varying  parameters  of  the  physical 
system . 


I V .  Simulation  Plan  C  3 ;  6 ;  93 

4.1.  Introduction 

Evaluation  of  the  performance  of  the  moving-bank 
multiple  model  adaptive  estimator/controller  for  this 
application  requires  simulating  actual  space  structure 
movement  and  est imator /control ler  operation.  The  computer 
simulations  provide  both  a  Monte  Carlo  analysis  and  a 
sensitivity  analysis  (using  ambiguity  functions)  of  the 
est imat or /cont r ol ler .  This  chapter  provides  background  on 
the  Monte  Carlo  simulation,  briefly  outlines  the  computer 
software,  and  then  discusses  the  simulation  plan  for 
analyzing  the  performance  of  the  estimator  design  logics 
and  the  moving-bank  algorithms. 

4.2.  Monte  Carlo  Ana  1  vs  Is 

The  objective  of  this  analysis  is  to  obtain 
statistical  information  on  the  est imator /control ler ' s 
performance,  through  the  use  of  a  Monte  Carlo  study.  A 
Monte  Carlo  simulation  must  be  used  Instead  of  a  covariance 
simulation,  even  though  the  system  is  linear,  because  of 
the  adaptive  nature  of  the  problem.  The  Monte  Carlo  study 
obtaining  many  samples  of  the  error  process  through 
simulation  and  then  using  this  data  to  approximate  ‘:he 
process  statistics  £21:3293. 

The  true  system  model  under  consideration  can  be 
described  by  a  linear  t Ime- i nvar lant  difference  equation: 


x(ti  +  i)  =  <|)(ti  +  i,ti  )x(ti  )  +  Bd(ti)u(ti)  +  Gd(ti)wd(ti) 

(4-1) 

(See  Equation  (2-1)  for  a  complete  definition  of  terms). 

and  are  the  discrete-time  equivalents  [22:17111  of 
the  B  and  G  matrices  given  in  Equation  (3-5).  It  is 
assumed  that  the  noise  input  matrix  is  identical  to  the 
control  input  matrix;  therefore  [21:1713, 

^d  =  I 

®d  ^  (t.  ,t)B  dr  (4-2) 

Noise-corrupted  measurements  are  provided  to  the  estimator 
in  the  form  of : 

z(t^)  =  Hx(tj^)  +  v(t^)  (4-3) 

where  H  is  the  measurement  matrix  and  v(t^)  is  a  discrete 
time,  zero-mean,  white  Gaussian  measurement  noise  with 
covariance  matrix  R.  Matrices  ♦,  B^,  G^,  and  H  are 
functions  of  the  true  parameter  vector  a^: 


where  and  are  the  true  mass  and  stiffness  para¬ 
meters,  respectively  and  are  discussed  in  Section  3.4. 


The  simulation  is  accomplished  for  a  sufficient 
number  of  runs  so  that  the  computed  sample  means  and 
variances  of  the  random  variables  of  interest  are  good 
approximations  to  true  ensemble  averages.  The  number  of 
simulation  runs  selected  is  10,  and  this  is  determined  as 
sufficient  by  observation  of  the  sample  statistics  as  the 
number  of  Monte  Carlo  runs  is  increased  C3:52J. 

Figure  4-1  illustrates  the  block  diagram  for  the 
Monte  Carlo  simulation  C3:53;  223.  The  variables  not 
previously  defined  are: 


x(t  ) 

i 

l(ti ) 


e  ( t .  ) 
^  1 


"truth  model"  states 
estimate  of  system  states 
estimate  of  uncertain  parameter  vector 
error  in  the  parameter  estimate 
e^Ct.  )  =  ^(t.  )  -  a(t.  ) 
error  in  the  system  state  estimate 


e  (t,  )  =  (t.  )  -  x(t,  ) 


where  T  is  a  n  x  n^  matrix  to  make  the  dimensions 
compatible,  since  the  estimate  is  typically  of  lesser 
dimension  than  that  of  the  "truth  model"  states. 


The  sample  mean  of  the  vector  of  variables  of 
Interest  Is  computed  as  C3:53;  6:46;  21:129;  9:613: 


System  Estimator  and  Controller  Simulation 


(4-4) 


ELe^(t,  )  ^ 
s:X  1 


M 

-ex 


(t,  )  =  (1/L)  E  (t  ) 

*  k=l  ^ 


where : 


L 


total  number  of  simulation  runs 

value  of  e_  (t.)  during  the  l<th  simulation  run 

X  i 


The  sample  covariance  of  e  (t.)  is  computed  as  C3:54; 

1 

6:47;  21:130;  9:631: 


L 

C1/(L-1)]  E 
k  =  l 


^k<‘l'Sxk“l' 


(4-5) 


Similar  relationships  for  e  (t.)  can  also  be  generated 
read i ly . 

When  evaluating  the  estimator  alone,  the  feedback 
controller  in  Figure  4-1  is  replaced  by  a  dither  signal 
with  a  frequency  and  amplitude  that  is  determined  by  trial 
and  error.  It  has  been  shown  that  a  dither  signal  can  be 
used  to  excite  the  system  model  and  enhance  parameter  iden 
tification  C6:50, 58-59;  20:135,136;  21:229;  91.  The 

error  in  the  state  estimate  and  the  error  in  the  parameter 
estimate  are  useful  in  evaluating  the  performance  of  the 
estimator.  The  error  in  the  state  estimate  gives 


the  best  means  of  comparing  the  estimator  to  other  types  of 
estimators,  since  the  primary  concern  is  accurate  state 
estimation  (and  control  eventually)  rather  than  parameter 
estimation  for  its  own  sa)<e.  The  error  in  the  parameter 
estimate  lends  insight  into  the  accuracy  of  state  esti¬ 
mation,  indicates  the  precision  of  a  parameter  estimate 
that  may  be  fed  to  an  adaptive  controller,  and  provides  a 
means  of  evaluating  various  move,  contract,  and  expand 
algorithms  C3:54]. 

When  evaluating  the  estimator/controller  combination, 
it  is  more  appropriate  to  examine  the  statistics  of  the 
true  state  values.  In  this  thesis  effort,  the  control 
objective  is  to  quell  any  oscillations  in  the  two-bay  truss 
structure  and  to  "point"  the  two-bay  truss  in  a  commanded 
direction.  It  is  also  important  to  examine  the  magnitude 
of  the  control  inputs  in  order  to  detect  unreasonable 
commanded  control  levels.  Lynch  I16D  in  his  research  with 
a  fixed  two-bay  truss  (see  Section  3.4.2)  limited  the 
magnitude  of  the  force  of  an  individual  actuator  to  100 
lbs.  Although  this  limit  is  not  implemented  in  this 
thesis,  it  is  used  as  a  general  guideline  as  to  a  reason¬ 
able  range  of  actuator  activity. 

4.3.  Software  Descr i pt i on  C3:55-593 

4.3.1.  Introduction.  The  analysis  of  the  moving 
-banVi  est  imator/controller  requires  the  development  of 


four  computer  programs.  Each  of  these  programs  is  a 
modification  of  programs  developed  by  Hentz,  Filios,  and 
Karnick  C3;  6;  91.  For  a  more  detailed  description  of  the 
following  programs,  the  reader  is  referred  to  Filios  C33. 
The  first  program  is  a  preprocessor  which  creates  a  para¬ 
meter  space  that  is  utilized  in  the  Monte  Carlo  simulation 
and  ambiguity  functions  analysis  (see  Section  2.5).  The 
second  program  simulates  the  moving-bank  multiple  model 
adaptive  estimator  and  performs  Monte  Carlo  simulation  runs 
and  generates  data  for  each  run.  This  program  is  also 
modified  to  produce  Monte  Carlo  runs  of  individual  ele¬ 
mental  filters.  This  allows  for  sensitivity  studies  that 
are  conducted  on  individual  controllers,  such  as  tuning  the 
state  weighting  matrix  (X)  and  control  weighting  matrix  (U) 
for  Individual  controllers.  The  third  program  is  a  post¬ 
processor  that  computes  the  means  and  variances  of  var¬ 
iables  of  interest  and  then  generates  the  plots  of  statis¬ 
tics  for  the  Monte  Carlo  simulation.  The  fourth  program 
computes  the  ambiguity  functions  and  generates  their  plots. 

4.3.2.  Preprocessor .  The  preprocessor  computes  the 
discretized,  reduced  order  (equal  to  the  truth  model), 
system  matrices  (  <*. ,  ,  G^,  H,  in  Equations  (2-1) 

and  (3-19))  for  each  parameter  point  within  the  para¬ 
meter  space,  the  Kalman  filter  and  LQ  controller  gains, 
and  information  needed  for  the  ambiguity  functions  analysis 


(see  Section  2.5).  An  input  file  allows  the  state  and 
control  weighting  matrices  (X  and  U),  dynamic  driving  noise 
strength  (Q),  measurement  noise  covariance  (R),  as  well  as 
the  time  increment  for  the  discrete  system,  to  be  varied. 
This  input  file  also  contains  the  mass  and  stiffness 
matrices  and  two  vectors  which  specify  the  mass  and  stiff¬ 
ness  variation.  The  mass  and  stiffness  matrices  are  used 
to  determine  a  24-state  system  (see  Equation  3-4)  which  is 
then  reduced  to  the  number  of  states  specified  (see  Section 
3.5). 

4.3.3.  Primary  Processor .  The  primary  processor 
performs  the  Monte  Carlo  simulations.  For  each  sample 
period,  the  true  system  and  the  filters  currently  imple¬ 
mented  in  the  moving  bank  are  propagated  forward  from  the 
most  recent  sample  time.  A  noise  corrupted  measurement  is 
then  made  of  the  true  system,  and  the  filters  of  the  moving 
bank  are  updated.  The  program  then  calculates  the  neces¬ 
sary  control  Inputs  and  makes  decisions  on  whether  to  move, 
expand,  or  contract  the  bank.  After  each  sample  period  is 
complete,  the  values  of  the  variables  of  interest  are 
written  to  a  data  file. 

The  inputs  to  the  primary  processor  describe  the 
parameter  space  (obtained  from  the  preprocessor)  and  true 
system  parameters,  and  specify  the  move/contract/expand 
algorithms  to  be  implemented,  the  associated  thresholds. 


primary  processor  is  a  data  file  for  each  variable  of 
interest  (state  estimate,  actual  state  values,  control 
inputs)  covering  all  of  the  simulation  runs,  and  a  more 
detailed  print  file  covering  just  the  first  Monte  Carlo 
run.  The  print  file  lists  the  exact  filters  implemented  in 
the  moving  bank  and  the  variables  which  affect  the  decision 
a  Igor  1 thms . 

4.3.4.  Postpr ocessor .  The  postprocessor  takes  the 
variable  data  files  obtained  from  the  primary  processor  and 
calculates  the  sample  means  and  variances  from  t^  to 
Plots  are  then  generated  which  contain  the  time  histories 
of  the  means  and  the  standard  deviation  of  each  variable. 
The  postprocessor  is  run  for  each  data  file  generated  by 
the  primary  processor. 


4.3.5. 


Functions  Ana  1 vs i s .  The  ambiguity 


function  analysis  involves  two  programs.  The  first 
evaluates  the  ambiguity  function  for  each  point  in  the 
parameter  space  (See  Chapter  II  for  a  full  discussion  of 
ambiguity  functions).  A  square  wave  dither  signal  is  used 
to  enhance  the  difference  between  filter  models.  The  form, 
magnitude,  and  frequency  of  the  dither  was  determined  by 
trial  and  error.  The  input  file  to  this  program  specifies 


the  truth  model,  number  of  Monte  Carlo  runs,  as  well  as  the 


time  length  of  the  simulation.  The  system  matrices  are 
obtained  from  the  same  file  generated  by  preprocessor  for 
the  primary  processor.  The  second  program  reads  the 
Information  from  a  file  and  then  generates  a  three- 
dimensional  plot  of  the  ambiguity  function  values  versus 
location  in  the  parameter  space. 

4.4.  Simulation  Plan 

4.4.1.  Noise  Level  Determination .  The  goal  of  the 
noise  level  determination  study  was  to  determine  an 
appropriate  value  for  the  noise  covariance  (R)  and  dynamics 
noise  strength  (Q).  Therefore,  the  Q  and  R  values  for  the 
"truth"  model  and  filter  model  were  set  equal  and  were 
assumed  to  be  diagonal  for  simplicity.  The  parameter  used 
for  both  the  "truth"  and  filter  model  was  (5,5),  which 
corresponds  to  the  center  of  the  bank.  This  particular 
parameter  was  chosen  because  it  is  the  center  of  the 
overall  parameter  space.  For  this  study,  the  single  filter 
primary  processor  consisted  of  a  single  filter  matched 
against  a  "truth"  model. 

The  Q  values  were  determined  by  holding  the  R  values 
constant  and  varying  the  Q  values  one  at  a  time.  The  Q 
values  were  decreased  by  orders  of  magnitude  until  the  rms 
error  for  the  state  estimate  no  longer  decreased  substan¬ 
tially.  After  determining  the  Q  values,  the  R  values  were 


found  using  the  same  technique  but  holding  the  Q  values 
constant  and  varying  the  R  values  one  at  a  time.  NoLe 
that,  during  this  study,  the  Q  and  R  values  for  the  "truth' 
and  filter  models  always  remain  equal. 

4.4.2.  Two-State  Model .  Hentz  showed  that  for  a 
two-state  (i.e.  one  bending  mode)  filter  model  matched 
against  a  two-state  truth  model,  the  center  of  the  moving 
bank  would  move  in  the  direction  of  the  true  parameter. 

The  goal  of  the  two-state  model  study  was  to  reproduce 
Hentz's  work  C61  for  the  truss  model.  To  accomplish  this, 
the  preprocessor  was  modified  so  that  it  produced  a  two- 
state  (i.e.,  one  bending  mode)  filter  for  each  parameter 
point.  The  primary  processor  in  the  two-state  model  study 
used  a  parameter  space  that  consisted  of  all  100  points 
upon  which  to  base  elemental  Kalman  filters.  In  addition 
the  primary  processor  was  modified  so  that  a  record  could 
be  kept  of  the  movement  of  the  estimated  true  parameter  (a) 
(i.e.  movement  of  the  bank),  throughout  the  parameter 
space.  The  record  was  acquired  by  initializing  a  10x10 
matrix  of  elements  to  zero  and  then  adding  unity  at  each 
sample  time  to  the  matrix  element  corresponding  to  the 
current  estimate  of  the  true  parameter.  Entries  were 
recorded  from  t=  1.0  to  the  end  of  the  run.  Recording  was 
not  started  at  t  =  0.0  in  an  attempt  to  avoid  transient 
effects.  In  addition,  the  moving-bank  used  the  smallest 


discretization  level  and  was  Initially  centered  at  the 
(5,5)  parameter  point  (i.e.,  the  middle  of  the  parameter 
space).  Finally,  the  postprocessor  was  not  used  in  this 
study  because  plots  of  the  state  variables  were  not 
required . 

4.4.3.  Six-State  Filter .  The  next  step  in  this 
research  effort  was  to  determine  if  system  complexity 
affects  the  performance  of  the  movlng-ban)c  algorithm.  The 
preprocessor  was  returned  to  its  original  form  of  producing 
six-state  filters.  The  filters  consisted  of  two  states  per 
mode  for  one  rigid  body  mode  and  two  bending  body  modes. 

The  primary  processor  required  no  modifications  from  the 
two-state  study.  Also,  the  postprocessor  was  not  needed 
for  this  study  because  the  plots  of  the  state  variables 
were  not  required. 

4.4.4.  Space  Discretization.  Following  the 
six-state  filter  study,  a  space  discretization  study  was 
conducted  to  determine  an  appropriate  space  discretization 
level.  The  preprocessor  again  produced  six-state  filters. 
However,  as  In  the  noise  level  determination  study,  the 
primary  processor  used  a  single  filter  in  the  banlt  matched 
against  a  "truth"  model.  Also,  the  postprocessor  was 
modified  to  produce  rms  error  plots.  The  space  discre¬ 
tization  study  was  conducted  using  the  following  steps; 


1.  choose  one  parameter  to  vary  and  hold  the  other 
parameter  constant  . 

2.  Set  the  filter  parameter  equal  to  the  middle  value 
in  the  range  of  possible  values  for  the  parameter  that  is 
constant  and  to  the  smallest  value  in  the  range  of  possible 
values  for  the  varying  parameter. 

3.  Set  the  truth  model  parameter  equal  to  the  filter 
parameter  value. 

4.  Obtain  RMS  error  plots  for  this  setup. 

5.  Set  the  truth  model  parameter  equal  to  the  next 
larger  parameter  value  for  the  varying  parameter. 

6.  Obtain  rms  error  plots  for  this  setup. 

7.  Compare  the  rms  error  plots  from  step  6  and  step 

4.  If  the  rms  error  plots  for  step  6  are  20%  to  30%  larger 

than  for  step  4,  leave  the  parameter  point  that  corresponds 
to  the  filter  in  the  parameter  space  and  move  the  filter  to 

the  parameter  point  associated  with  the  truth  model  (repeat 

steps  2  through  7).  If  the  rms  error  plots  for  step  6  are 
less  than  20%  larger  than  the  plots  in  step  4,  remove  the 
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parameter  point  associated  with  the  truth  model  from  the 
parameter  space  and  move  the  truth  model  to  the  next  larger 
parameter  in  parameter  space  (repeat  steps  5  through  7). 

If  the  rms  error  plots  for  step  6  are  greater  than  30% 
larger  than  the  plots  In  step  4,  move  the  truth  model 
halfway  between  the  filter  parameter  and  the  truth  model 
parameter  (repeat  steps  5  through  7). 


8.  Interchange  the  varying  parameter  and  the  constant 
parameter  after  the  varying  parameter  space  is  determined. 
Then  repeat  the  process  for  the  other  parameter. 


Note:  If  the  edge  of  parameter  space  Is  reached 

before  10  values  for  the  parameter  are  determined,  then 
augment  the  parameter  space  with  larger  parameter  values 


The  objective  of  this  study  was  to  obtain  an  rms  error 
growth  of  20%  to  30%  between  each  step  In  parameter  space 
The  results  of  this  study  are  contained  in  Section  5.5, 
with  the  eventual  ranges  of  the  mass  and  stiffness 
parameters  shown  In  Figure  5-11. 

4.4.5.  Measurement  and  Dynamics  Noise 
Sensitivity.  The  goal  of  this  study  was  to  determine  the 


upper  limits,  whose  existence  Is  discussed  in  Section  1.4, 
on  the  Q  and  R  values.  Once  again,  the  prepr occessor 


produced  six-state  filters  and  the  primary  processor  used  a 
single  filter  matched  against  a  truth  model.  As  in  the 
noise  level  determination,  the  measurement  noise  covari¬ 
ances  (R)  and  the  dynamics  noise  strengths  (Q)  for  the 
"truth"  model  and  the  filter  model  were  set  equal  to  each 
other.  The  Q  limits  were  determined  by  holding  the  R 
values  constant  while  increasing  the  Q  values  one  at  a 
time.  Each  Q  value  was  increased  by  orders  of  magnitude 
until  the  moving-bank's  estimate  no  longer  locked  onto  the 
true  parameter.  After  the  Q  limits  were  found,  the  R 
limits  were  determined  by  holding  the  Q  values  constant  and 
increasing  the  R  values,  one  at  a  time,  by  orders  of 
magnitude  until  the  moving-bank's  estimate  of  the  true 
parameter  no  longer  locked  onto  the  actual  true  parameter. 
This  degradation  in  performance  was  determined  by  main¬ 
taining  a  record  of  the  true  parameter  estimate  as  in  the 
two-state  filter. 

4.4.6.  Controller  Design ;  Need  for 

Cross  Terms  in  the  Quadratic  Cost 

The  goal  of  this  study  was  to  determine  if 
the  moving-bank  controller  performance  could  be  increased 
by  adding  cross  terms  to  the  quadratic  cost  function.  The 
cross  weighting  matrix  (S)  was  evaluated  L23:763  and  com¬ 
pared  to  the  state  weighting  matrix  (X)  and  the  control 
weighting  matrix  (U).  The  generalized  cost  function  used 
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in  this  study  is  given  in  Equation  (4-6)  t23:733.  The 
differential  equations  used  to  solve  for  the  cross  terms 
are  given  in  (4-7)  through  (4-12). 


J  =  E|  (l/2)x'^(tj^^^)XjX(tj^^j^)  +  £(l/2)Cx'^(t.  )X(t^  )x(t^  ) 


+  ^"^(t.  )U(t.  )u(tp  +  2x'^(tj  )S(t.  )u(tpif 


(4-6) 


where 


X  ( t^  ) 

u(t.  ) 
X(t^  ) 
U(t.  ) 
S(t.  ) 


=  state  vector 


=  control  input  vector 
=  state  weighting  matrix 
=  control  weighting  matrix 


=  cross  term  matrix 


♦  (t,t.  ) 
B(t,t.  ) 
X(t,t.  ) 


=  F(t)  <|>(t,t^) 

=  F(t)i(t,t. )  +  B(t) 

=  <l>'^(t,t^  ^ 


Q(t,t.)  =  F(t)Q(t,t.)  +Q(t,t. )F^(t)  + 


G(t)Q(t)GMt) 


(4-7) 

(4-8) 

(4-9) 


(4-10) 


S(t,t.  ) 


U(t,t, 


=  4.'^(t,t^  )W^^(t)B(t,t.  )  +  ♦'^(t^t^  )W^^_|(t) 

(4-11) 


=  b' (t,t.  )W^^(t)B(t,t.  )  +  Wj^^(t) 

+  B’’(t,tj)W^^(t)  +  wJ^^(t)B(t,t^) 


(4-12) 


The  state  transition  matrix,  ♦  ,  is  initially  set  equal  to 

the  identity  matrix,  while  all  the  other  initial  conditions 

are  set  equal  to  a  properly  dimensioned  zero  matrix.  Also 

the  W  matrix  was  assumed  to  be  equal  to  0.  These  dif- 
xu 

ferentlal  equations  were  solved  using  MACSYMA.  In  addi¬ 
tion,  because  the  matrices  were  of  different  dimensions, 

T 

SS  was  formed  and  then  the  square  root  of  the  eigen- 

T 

values  of  the  SS  matrix  were  used  to  provide  a  compari- 
parlson  with  the  eigenvalues  of  the  X  matrix.  During  this 
study,  it  was  determined  that  the  cross  terms  do  not  sig¬ 
nificantly  contribute  to  the  generalized  cost  function. 

The  detailed  results  of  this  study  are  contained  in  Section 


4.4.7.  State  Weighting  and  Control  Weighting 
Matrix  Tuning .  The  object  of  this  study  was  to  determine 
appropriate  values  for  the  X  and  U  matrices.  The 
prepr occessor  again  produced  six-state  filters,  while  the 
primary  processor  used  a  single  filter  bank  and  was  modi¬ 
fied  to  output  the  true  state  vector.  The  true  parameter 
and  the  filter  parameter  were  both  set  to  (7,6).  The 
parameter  point  (7,6)  was  chosen  so  that  it  was  near  the 
center  of  the  bank  but  not  equal  to  a  parameter  point 
corresponding  to  one  of  the  fixed-bank  filters  C(l,l), 
(1,5),(1,9),(5,1),(5,5),(5,9),(9,1),(9,5),(9,9)J.  The 
state  weighting  matrix  (X)  values  were  determined  by 
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holding  the  control  weighting  matrix  (U)  constant  and 
increasing  the  X  values  one  at  a  time  until  the  rms  values 
for  the  true  states  stop  decreasing  drastically  as  the  X 
values  were  increased.  After  the  X  values  were  determined, 
the  U  values  were  found  by  holding  the  X  values  constant 
and  decreasing  the  U  values  until  the  rms  values  for  the 
true  states  stop  decreasing  drastically  as  the  U  values 
were  decreased. 

4.4.8.  Controller  Study.  The  goal  of  the  controller 
study  was  to  compare  the  performance  of  a  moving-bank 
controller  with  a  fixed-bank  controller  and  a  benchmark 
controller.  The  prepr occessor  again  produced  six-state 
filters,  while  the  parameter  space  used  in  the  primary 
processor  consisted  of  the  entire  100  points  and  the 
primary  processor  was  modified  to  produce  the  true  state 
vector .  The  postprocessor  produced  the  true  state  vector 
and  control  vector  statistics  plots.  This  setup  was  used 
for  the  moving-bank  controller  and  the  fixed-bank  con¬ 
troller.  However  in  producing  the  benchmark  controller, 
the  filter  parameter  and  true  parameter  values  were  set 
equal  to  each  other  and  the  parameter  space  used  in  the 
primary  processor  consists  of  one  (true)  point  upon  which 
to  base  one  Kalman  filter.  Thus  the  optimal  control  was 
produced  because  the  filter  knows  the  true  parameter.  The 
comparison  deals  mainly  with  each  controller's  ability  to 
quell  oscillations. 
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4.5.  Summary 

The  analysis  technique  used  in  this  thesis  is 
described  in  Section  4.2.  In  Section  4.3,  the  simulations 
themselves  are  described.  Finally,  the  simulation  plan 
used  in  this  thesis  is  detailed  in  Section  4.4. 
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V.  Results 


5.1.  Introduction 

This  chapter  presents  the  results  of  the  Monte  Carlo 
simulations.  The  goal  of  this  research  effort  was  to 
ascertain  the  cause  behind  the  poor  performance  of  the 
moving  bank  in  Karnlck's  research  C93  and  then  to  determine 
the  optimal  moving-bank  or  fixed-bank  structure  in  order  to 
provide  the  best  control  of  the  large  space  structure 
model.  The  first  step  was  to  find  the  dynamics  noise 
strength  (Q)  and  the  measurement  noise  covariance  (R)  that 
provided  for  low  rms  error  on  the  state  variables,  while 
not  allowing  the  Q  and  R  values  to  decrease  to  zero:  a 
point  of  diminishing  returns  was  subjectively  established. 
Recall  from  Section  4.4.1,  that  this  was  accomplished  by 
setting  the  "truth"  and  filter  Q  and  R  values  equal  to  each 
other.  These  Q  and  R  values  were  then  varied  one  at  a 
time.  With  these  appropriate  Q  and  R  values,  the  center  of 
the  moving  bank  was  pulled  in  the  direction  of  the  true 
parameter,  rather  than  wandering  aimlessly  in  parameter 
space.  The  system  complexity  was  then  Increased  to  six 
states  to  determine  if  the  performance  of  the  moving  bank 
decreased  as  the  system  complexity  increased.  From  this 
study,  it  was  determined  that  the  performance  of  the  moving 
bank  did  not  degrade  as  system  complexity  Increased. 
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Therefore,  the  problem  behind  the  poor  performance  of  the 
moving  bank  was  the  choice  of  the  Q  and  R  values:  there  Is 
a  fundamental  requirement  on  the  range  of  the  Q  and  R 
values  above  which  the  adaptive  mechanization  of  the 
moving  bank  breaks  down.  Section  5.6  discusses  the  range 
requirements  on  the  Q  and  R  values  for  this  model. 

However,  even  though  the  moving  bank  supplied  an 
accurate  estimate  of  the  true  parameter,  some  of  the  state 
estimation  error  plots  were  unstable.  In  an  attempt  to 
correct  this  problem  a  space  discretization  study  was 
performed.  After  the  completion  of  the  space  discreti¬ 
zation  study,  the  state  estimation  error  plots  were  then 
stable  (see  Section  5.5).  Subsequently,  a  study  was 
conducted  to  determine  if  additional  performance  could  be 
attained  by  adding  cross  terms,  between  states  and  con¬ 
trols,  to  the  quadratic  cost  function  used  in  the  con¬ 
troller  synthesis.  Section  5.7  shows  that  the  addition  of 
cross  terms  to  the  quadratic  cost  function  did  not  appre¬ 
ciably  enhance  the  controller  performance.  After  consid¬ 
ering  cross  terms,  additional  performance  was  achieved  by 
tuning  the  state  weighting  matrix  (X)  and  the  control 
weighting  matrix  (U).  This  tuning  was  performed  for  the 
(7,6)  parameter  point,  which  was  chosen  so  that  the  tuned 
parameter  point  would  be  near  the  center  of  the  bank 


but  not  equal  to  a  parameter  value  corresponding  to  one  of 
the  fixed  bank  filters  C ( 1, 1 ) , ( 1, 5 ) , ( 1, 9 ) , ( 5, 1 ) , ( 5, 5 ) , ( 5, 

9  ) , ( 9, 1 )  , ( 9, 5  ) , ( 9, 9  )J  .  The  results  of  this  tuning  process 
are  given  in  Section  5.8.  The  final  attempt  to  enhance 
performance  was  the  use  of  a  transition  period  of  partially 
adaptive  control  between  the  time  used  to  estimate  states 
in  open  loop  fashion  (with  a  dither  signal  to  excite  the 
system)  and  the  time  when  full  multiple  model  adaptive 
control  was  applied  to  the  system.  It  was  thought  the 
transition  period  was  needed  because  if  adaptive  control 
inputs  were  supplied  to  the  system  before  an  accurate 
estimate  of  the  true  parameter  was  obtained,  the  control 
inputs  could  drive  the  system  unstable.  Therefore,  a 
nonadaptive  controller  was  designed  that  was  very  robust  to 
parameter  variation  at  the  expense  of  controller  perform¬ 
ance  so  that  the  controller  inputs  would  at  least  maintain 
stability  regardless  of  the  estimate  of  the  true  parameter. 
However,  Section  5.9  shows  that  this  transition  period  was 
not  needed  because  the  moving-bank  algorithm  locks  onto  the 
true  parameter  in  less  than  0.5  sec,  thereby  allowing  full 
adaptive  control  to  be  applied  to  the  system  immediately 
following  the  excitation  and  open  loop  estimation  period. 

It  was  also  determined  during  this  controller  evaluation, 
the  state  weighting  matrix  (X)  and  the  control  weighting 
matrix  (U)  needed  to  be  tuned  for  each  parameter  point 


in  order  to  achieve  the  best  performance.  This  is  also 
covered  in  Section  5.9. 

Finally  after  all  the  techinques  used  to  enhance 
performance  were  applied  to  the  system,  a  performance 
comparison  between  a  benchmark  controller,  a  moving-bank 
controller,  and  a  fixed-bank  controller  was  conducted,  with 
the  results  shown  in  Section  5.9.  This  performance  com¬ 
parison  showed  that  the  moving-bank  algorithm  converges  to 
the  performance  of  the  benchmark  controller,  while  the 
fixed-bank  controller  converges  to  the  benchmark  perform¬ 
ance  only  when  the  true  parameter  corresponds  to  one  of  the 
filters  Implemented  in  the  fixed-bank  controller.  The 
benchmark  controller  is  described  in  Section  5.9  and  Is 
considered  to  provide  the  optimal  control  performance. 

5.2.  Measurement  Covariance  Noise  and  Dynamics  Noise 

Strength  Determination 

Before  the  research  can  be  started,  appropriate 
choices  of  the  measurement  noise  covariance  (R)  and 
dynamics  noise  strength  (Q)  must  be  found.  A  description 
of  the  process  used  to  determine  appropriate  choices  of  the 
R  and  Q  values  is  given  in  Section  4.4.1.  Recall  though 
that  the  truth  and  filter-assumed  values  were  set  equal  to 
each  other,  then  both  the  truth  and  filter  values  for  the 
entries  in  these  matrices  are  varied  one  at  a  time.  Also 
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note,  that  the  measurement  noise  covariance  (R)  and 
dynamics  noise  strength  (Q)  are  determined  using  a  single 
filter  primary  processor.  The  resulting  Q  and  R  values 
obtained  from  this  process  are  given  in  Table  5-1.  Table 
5-2  contains  the  time  averaged  state  estimation  error 
statistics,  omitting  the  initial  transient  time,  for  the  Q 
and  R  values  given  in  Table  5-1.  These  statistics  were 
also  produced  using  Q  and  R  values  which  were  one  order  of 
magnitude  larger  than  in  Table  5-1.  These  time  averaged 
statistics  are  contained  in  Table  5-3.  As  can  be  seen,  the 
rms  errors  in  Table  5-3  are  about  one  order  of  magnitude 
larger  than  those  in  Table  5-2.  Therefore,  the  Q  and  R 
values  contained  in  Table  5-1  provide  better  estimation 
performance  than  the  Q  and  R  values  used  in  Table  5-3. 
After  appropriate  Q  and  R  values  are  determined,  the 
remainder  of  the  research  effort  can  be  performed. 

5.3.  Two-State  Mode  1 

The  goal  of  the  two-state  model  was  to  determine  if 
Hentz's  work  C63  could  be  duplicated  for  the  large  space- 
structure  model.  Hentz  showed  that,  for  a  two-state  (i.e. 
one  bending  mode)  filter  model  matched  against  a  two-state 
truth  model,  the  center  of  the  moving  bank  would  move  in 
the  direction  of  the  true  parameter.  However,  it  is 
important  to  note  that  optimal  performance  of  the 


Table  5~1 


Dynamics  and  Measurement  Noise  Strength 
(for  both  the  truth  and  filter  models) 

Dynamics  Noise  Strengths 

2  2  3 

(Units:  Position  -  in  /sec;  Velocity  -  in  /sec  ) 

Noise  Input  at  Node  1  (Ql)  -  80. 

Noise  Input  at  Node  2  (Q2)  -  80. 

Noise  Input  at  Hub  (Q3)  -  50. 

at  angular  velocity 
level 

Noise  Input  at  Hub  (Q4)  -  50. 
at  angular 
acceleration  level 


Measurement  Noise  Variances 
2  2 

(Units:  Position  -  in  sec;  Veloctiy  -  in  /sec) 


Pos i t i on 

Sensor 

at 

Node 

1 

(Rl) 

-  0.000027 

Ve loc i ty 

Sensor 

at 

Node 

1 

(R3) 

-  0.0025 

Position 

Sensor 

at 

Node 

2 

(R2) 

-  0.0000027 

Velocity 

Sensor 

at 

Node 

2 

(R4) 

-  0.0042668 

Position 

Sensor 

at 

Hub 

(R5) 

-  12.0 

Velocity 

Sensor 

at 

Hub 

(R6) 

-  0.000855 

Table  5-2 

Time  Averaged  State  Estimation  Errors  Statistics: 
Using  the  Q  and  R  Values  Given  in  Table  5-1; 

For  Parameter  Point  (5,5) 

(Units:  Position  -  in;  Velocity  -  In/sec) 


State 


Standard  Deviation 


Position 

Error 

(Node 

1) 

0.008639 

0.002310 

Velocity 

Error 

( Node 

1) 

0.05540 

0.006272 

Pos i t i on 

Error 

(Node 

2) 

0.01925 

0 . 0008319 

Velocity 

Error 

( Node 

2) 

0.03963 

0 .046079 

Pos i 1 1  on 

Error 

( Hub ) 

-0.21341 

0.43714 

Velocity 

Error 

(Hub) 

0.9989 

3.3993 

Table  5-3 

Time  Averaged  State  Estimation  Errors  Statistics: 
Using  the  Q  and  R  Values  One  Order  of  Magnitude  Larger 

than  in  Table  5-1; 

For  Parameter  Point  (5,5) 

(Units:  Position  -  in;  Velocity  -  in/sec) 


State 


Standard  Deviation 


Pos i t ion 

Error 

(Node 

1) 

0.0202 

0.023 

Velocity 

Error 

( Node 

1) 

0.172 

0.06284 

Position 

Error 

( Node 

2) 

0.0657 

0 . 00822 

Velocity 

Error 

(Node 

2) 

0.12559 

0.4607 

Position 

Error 

( Hub ) 

-0.673 

4 . 3713 

Velocity 

Error 

(Hub) 

3.1657 

33.993 

moving  bank  was  not  a  goal  of  this  study.  Recall  from 
Section  4.4.2,  that  both  the  truth  and  filter  models  were 
reduced  to  two  states  (i.e.,  one  bending  mode).  In  the 
Monte  Carlo  simulations,  a  dither  signal  (square  wave)  with 
a  magnitude  of  300  and  frequency  of  30  rad/sec  was  used  to 
excite  the  system.  The  magnitude  and  frequency  of  the 
dither  signal  were  chosen  by  trial  and  error  using  the 
crltiera  that  the  system  be  excited  but  not  dominated  by 
the  dither  signal.  In  order  to  determine  the  trends  in  the 
movement  of  the  bank,  a  record  was  needed  of  the  estimated 
true  parameter  C93.  This  record  was  acquired  by  initial¬ 
izing  a  10x10  matrix  of  elements  to  zero  and  then  adding 
one  at  each  sample  time  to  the  matrix  element  corresponding 
to  the  estimate  of  the  true  parameter.  Entries  were 
recorded  from  t  =  1,0  to  the  end  of  the  run.  Recording  was 
not  started  at  t  =  0.0  in  an  attempt  to  avoid  transient 
effects.  The  transient  effects  could  be  omitted  because 
they  were  short  lived  and  stable.  This  is  shown  in  Section 
5.9.  In  addition,  the  moving-bank  algorithm  used  the 
smallest  discretization  level  and  was  initially  centered  at 
the  (5,5)  parameter  point  (i.e.,  in  the  middle  of  the  para¬ 
meter  space).  Using  the  above  techniques,  data  for  Figure 
5-1  through  Figure  5-3  was  acquired  for  three  true  para¬ 
meter  points.  In  Figure  5-1,  the  estimated  true  parameter 
was  spread  out  somewhat  in  parameter  space.  However,  most 
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of  the  data  points  were  contained  In  the  lower  left  quad¬ 
rant  which  indicates  a  pulling  of  the  center  of  the  moving 
bank  in  the  correct  direction.  In  Figure  5-2,  the  mass 
parameter  is  pulled  in  the  correct  direction  but  the 
stiffness  parameter  is  not  pulled  in  any  direction.  Again 
in  Figure  5-3,  the  mass  parameter  is  pulled  in  the  correct 
direction  while  this  time  the  stiffness  parameter  is  actu¬ 
ally  pulled  in  the  wrong  direction.  However,  there  is  an 
overall  indication  of  a  pulling  trend  of  the  center  of  the 
moving  bank  toward  the  true  parameter.  Therefore,  Hentz's 
work  was  duplicated  for  this  model.  The  quality  of  the 
state  estimate  produced  by  the  moving-bank  algorithm  is 
given  in  the  six-state  filter  study. 

5.4.  Six-State  Filter 

Since  Hentz's  work  was  duplicated  for  the  one-bending 
mode  model,  it  was  desired  to  determine  if  the  performance 
of  the  moving-bank  algorithm  degrades  for  a  more  complex 
system.  Therefore,  the  state  dimension  of  both  the  filter 
model  and  the  "truth"  model  were  increased  to  six  (i.e., 
one  rigid  body  mode  and  two  bending  body  modes).  Recall 
from  Section  4.4.3,  there  was  no  purposeful  order  reduction 
in  generating  the  filter  model,  and  thus  the  robustness  of 
the  algorithm  to  such  model  mismatch  has  not  as  yet  been 
explored.  The  data  shown  in  Figures  5-4  through  5-7  was 


acquired  using  the  same  techniques  as  the  data  used  In 
Figures  5-1  In  the  two-state  model  study  given  In  Section 
5.3,  except  the  model  now  has  six  states.  Figures  5-4 
through  5-7  show  that  there  Is  a  definite  pulling  of  the 
center  of  the  moving  bank  toward  the  true  parameter  point. 
In  analyzing  Figures  5-4  through  5-7  as  compared  to  Figures 
5-1  through  5-3,  It  appears  that  the  six-state  system 
performed  better  than  the  two-state  system,  which  would 
Indicate  that  the  moving  bank  performed  better  as  the 
system  complexity  Increases.  However,  recall  from  Section 
5.2,  that  the  Q  and  R  values  were  chosen  for  the  six-state 
model.  These  Q  and  R  values  were  used  In  the  two-state 
model  study,  therefore  the  Q  and  R  values  were  not  tuned 
fir  the  two-state  model  study.  However,  the  goal  of  the 
two-state  model  study  was  to  show  the  trend  of  the  center 
of  the  bank,  not  to  obtain  optimal  performance.  Therefore 
the  six-state  model  should  outperform  the  two-state  model. 
The  performance  of  the  six-state  filter  Indicates  that 
system  complexity  does  not  appreciably  affect  the  perform¬ 
ance  of  the  moving-bank  parameter  estimation.  However,  the 
six-state  model  still  has  a  performance  problem.  This 
problem  Is  Illustrated  In  Figures  5-8  through  5-10,  which 
contain  the  state  estimation  error  for  the  (2,7),  (9,9), 
(7,6)  parameter  points.  In  these  figures,  the  mean  value 
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(symbol  x),  the  mean  plus  one  sigma  (symbol  +),  and  the 
mean  minus  one  sigma  (symbol  y)  are  plotted  for  each  of  the 
six  states  in  the  filter.  Recall  In  Figure  5-4,  the 
moving-bank's  estimate  of  the  true  parameter  did  not  lock 
onto  the  actual  true  parameter  value.  The  corresponding 
state  estimation  errors  for  the  (2,7)  parameter  point  are 
unstable  (see  Figure  5-8).  However,  recall  In  Figure  5-6, 
the  moving  bank's  estimate  of  the  true  parameter  did  lock 
onto  the  actual  true  parameter  value.  The  corresponding 
state  estimate  error  plots  for  the  (9,9)  parameter  point 
are  stable  (see  Figure  5-9).  This  seems  to  Indicate  that 
there  is  a  strong  correlation  between  the  ability  of  the 
moving-bank  algorithm  to  estimate  the  true  parameter  and 
its  ability  to  provide  precise  state  estimates.  Couple 
this  correlation  with  the  fact  that  the  old  space  discre¬ 
tization  level  (see  Figure  5-11)  was  determined  by  intu¬ 
ition.  This  indicates  that  a  space  discretization  level 
study  must  be  performed  to  determine  an  appropriate  space 
discretization  level. 


5.5.  Space  Discretization 

Recall  from  Section  4.4.4,  the  space  discretization 
study  was  accomplished  by  monitoring  the  rms  error  of  the 
state  variables  as  the  true  parameter  and  the  filter  para¬ 
meter  are  moved  apart.  The  results  of  this  study  are 


given  in  Figure  5-11.  Also  note,  that  both  the  Inconsist¬ 
ent  spacing  and  the  enlarged  range  of  the  stiffness  and 
mass  values  were  chosen  to  enhance  parameter  identlfl- 
ablity.  After  the  optimal  space  discretization  level  was 
chosen.  Figure  5-12  through  Figure  5-15  were  produced  to 
compare  with  Figure  5-4  through  Figure  5-7.  Again  these 
figures  were  produced  using  the  same  technique  as  the 
figures  in  the  two-state  model  study  in  Section  5.3  except 
this  model  contains  six  states.  The  benefits  of  the  new 
space  discretization  level  can  be  seen,  especially  in 
Figure  5-12  when  compared  to  Figure  5-4.  In  Figure  5-12, 
the  moving-bank  estimate  of  the  true  parameter  locked  onto 
the  actual  true  parameter,  while  in  Figure  5-4,  the 
moving-bank  estimate  of  the  true  parameter  did  not  lock 
onto  the  actual  true  parameter.  In  addition,  plots 
containing  the  state  estimation  errors  were  produced  for 
the  true  parameters  (2,7),  (9,9),  and  (7,6)  for  the  new 
space  discretization  level  (see  Figures  5-16  through  5-16). 
In  comparing  Figure  5-8  with  Figure  5-16,  the  moving-bank 
algorithm  estimate  of  the  true  parameter  locked  onto  the 
actual  true  parameter  (see  Figure  5-12)  and  the  error  in 
the  state  estimate  Indicates  the  system  is  stable  (see 
Figure  5-16),  while  In  Figure  5-0  the  moving-bank  algorithm 
estimate  of  the  true  parameter  did  not  lock  onto  the  actual 


true  parameter  and  the  error  in  state  estimate  indicates 
that  the  system  is  unstable. 
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Fig.  5-12. 
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5.6.  Measurement  and  Dynamics  Noise  Sensitivit 


Since  the  choice  of  the  measurement  noise  covariance 
(R)  and  dynamics  noise  strength  (Q)  had  such  a  dramatic 
effect  on  the  performance  of  the  moving  bank  in  parameter 
estimation  (i.e.,  if  the  Q  and  R  values  were  chosen  too 
large  the  moving-bank  algorithm  fails  to  operate  cor¬ 
rectly),  a  study  was  performed  to  determine  when  the 
moving-bank  adaptive  controller  performance  degraded  to  an 
unacceptable  level  as  a  function  of  the  measurement  noise 
covariance  and  dynamics  noise  strengths.  Recall  from 
Section  4.4.5,  that  the  study  was  conducted  by  setting  the 
Q  and  R  values  for  the  truth  and  filter  models  equal  and 
then  increasing  these  Q  and  R  values  one  at  a  time.  The 
results  of  this  study  are  contained  in  Table  5-4.  For  some 
values  of  Q  and  R,  the  moving-bank  algorithm  showed  no  loss 
of  adaptive  performance  even  after  the  values  were 
increased  by  eight  orders  of  magnitude.  These  values  are 
listed  in  Table  5-4  as  showing  no  loss  of  adaptive  per¬ 
formance.  As  is  shown  in  the  table,  only  three  values 
significantly  affected  the  performance  of  the  moving  bank 
in  parameter,  and  these  values  are  all  contained  in  the 
measurement  covariance  matrix.  However,  it  makes  sense 
that  the  moving  bank  would  be  more  sensitive  to  R  values 
because  noisy  measurement  values  would  tend  to  mask  the 
difference  between  good  and  bad  assumed  models  within 


Table  5-4 


Dynamics  Noise  Strength  and  Measurement  Noise  Covariance 
Maximum  Values  for  Acceptable  Moving-Bank  Performance 


Dynamics  Noise  Strength 

2  2  3 

(Units:  Position  -  in  /sec;  Velocity  -  in  /sec  ) 

Noise  Input  at  node  1  (Ql)  -  No  loss  of  adaptive 

performance 

Noise  Input  at  node  2  (Q2)  -  No  loss  of  adaptive 

per f ormance 


Noise  Input  at  Hub  (Q3)  -  No  loss  of  adaptive 
at  angular  velocity  performance 

leve  1 


Noise  Input  at  Hub  (Q4)  -  No  loss  of  adaptive 
at  angular  acceleration  performance 

leve  1 


Measurement  Covariance  Noise 

2  2 

(Units:  Position  -  in  sec;  Velocity  -  in  /sec) 


Pos i t i on 

Sensor 

at 

Node 

1 

(Rl) 

-  0.27 

Pos i t i on 

Sensor 

at 

Node 

2 

(R2) 

-  0.027 

Velocity 

Sensor 

at 

Node 

1 

(R3) 

-  No  loss 

of 

adaptive 

per  f ormance 

Velocity 

Sensor 

at 

node 

2 

(R4) 

-  4.2668 

Pos i t i on 

Sensor 

at 

Hub 

(R5) 

-  No  loss 

of 

adaptive 

performance 


Velocity  sensor  at  Hub 


(R6)  -  No  loss  of  adaptive 
performance 


elemental  filters,  making  It  hard  to  discern  the  effect  of 
the  system  parameters.  Figure  5-19  through  Figure  5-22 
were  produced  using  the  same  technique  as  in  the  figures  in 
Section  5.3  with  Rl,  R2,  and  R4  simultaneously  set  to  the 
values  given  in  Table  5-4.  These  figures  show  that  the 
estimate  of  the  true  parameter  was  beginning  to  degrade,  as 
compared  to  the  performance  seen  in  Figures  5-12  through 
5-15.  For  example,  in  Figure  5-19  the  moving-bank  estiotate 
of  the  true  parameter  does  not  lock  onto  the  actual  true 
parameter,  while  in  Figure  5-12  the  estimate  of  the  true 
parameter  does  lock  onto  the  actual  true  parameter.  To 
show  that  this  degradation  continues  to  the  point  were  the 
moving  bank  completely  fails  to  operate.  Figure  5-23  and 
Figure  5-24  were  generated.  Again  these  figures  were 
generated  using  the  same  techniques  as  in  figures  in 
Section  5.3.  In  both  Figure  5-23  and  Figure  5-24,  Rl,  R2, 
and  R4  were  increased  an  additional  3  orders  of  magnitude 
beyond  the  values  in  Table  5-4.  Both  figures  show  that  the 
moving  bank  does  not  move  off  the  (5,5)  parameter  point, 
which  is  where  the  moving  bank  is  initially  centered. 


Therefore,  there  were  upper  limits  on  the  R  values,  for 
viable  parameter  estimation.  These  upper  limits  are 
Important  because  they  determine  the  precision  level  the 
sensors  must  meet  for  the  moving-bank  algorithm  to  operate 
correctly. 
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-state  Model  Bank  Location  Time  History 
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Fig.  5-21.  Six-State  Model  Bank  Location  Time  History 
True  Parameter  at  Mass  =  9,  Stiffness  =  9 
R1,R2,  and  R4  Set  to  Maximum  Value 


Fig.  5-22.  Six-State  Model  Bank  Location  Time  History 
True  Parameter  at  Mass  =  3,  Stiffness  =  3 
Rl,R2,  and  R4  Set  to  Maximum  Value 
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Rl,R2,  and  R4:  Three  Orders  of  Magnitude 
Larger  Than  In  Table  5-2 


Six-State  Model  Bank  Location  Time  History; 
True  Parameter  at  Mass  =  9,  Stiffness  =  9 
R1,R2,  and  R4:  Three  Orders  of  Magnitude 
Larger  Than  In  Table  5-2 
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Controller  Design :  Need  for 


Cross  Terms  In  the  Quadratic  Cost 
After  the  MMAE  began  to  supply  accurate  state 
estimation,  a  study  was  conducted  to  determine  if  addi¬ 
tional  performance  could  be  attained  by  adding  cross  terms 
to  the  quadratic  cost  function.  The  generalized  cost 
function  Is  given  in  Equation  (5-1)  C23:733.  The  LQ 
controller  Is  developed  fully  in  Appendix  A. 


J 


E  (l/2)x'^{t^^j^)X^x(tj^^j^)  +  £(l/2)Cx’^{tj  )X(tj  )x(t^  ) 

+  ^"^(t.  )U(tpy(t^  )  +  2x'^(tj^  )S(t^  )u(t^  )  ]  (5-1) 


where 


x(t^)  =  state  vector 
u(t^)  =  control  input  vector 
X(tj^)  =  state  weighting  matrix 
Utt^)  =  control  weighting  matrix 
S(tj^)  =  cross  term  matrix 


I 

i 

‘  The  cross  term  arises  due  to  the  fact  that  control  is  being 

applied  over  the  entire  sample  period,  not  just  at  the 
discrete  sample  tiroes  C23;75i.  Recall  that  the  differ¬ 
ential  equations  given  in  Section  4.4.6  C22:76D  were 

solved  to  determine  the  state  weighting  matrix  (X),  the 
control  weighting  matrix  (U),  and  the  cross  term  weighting 


matrix  (S).  In  addition  because  the  matrices  were  of  dif- 

T 

ferent  dimensions,  SS  was  formed  and  then  the  square 

T 

root  of  the  eigenvalues  of  the  SS  matrix  was  used  to  pro¬ 
vide  a  comparison  with  the  eigenvalues  of  the  X  matrix. 

T 

The  eigenvalues  for  the  U,  X,  and,  SS  matrices  are  given 

in  Table  5-5.  This  table  shows  that  the  eigenvalues  of 
T 

the  SS  matrix  are  at  least  an  order  of  magnitude  smaller 
than  the  eigenvalues  of  the  X  matrix.  Therefore,  the  cross 
weighting  matrix  does  not  attribute  greatly  to  the  cost 
function  and  was  not  added  to  the  quadratic  cost  function 
for  the  remainder  of  the  research.  However,  if  the  sample 
period  were  chosen  smaller,  then  the  cross  weighting  matrix 
would  have  a  larger  effect  on  the  cost  function.  The 
corresponding  state  control  performance  will  be  seen  in 
Section  5.9. 


5.8.  State  Weighting  and  Control  Weighting  Matrix  Tuning 
After  the  decision  was  made  not  to  include  a  cross 
weighting  matrix,  a  study  was  performed  to  tune  the  state 
weighting  matrix  (X)  and  the  control  weighting  matrix  (U). 
Recall  from  Section  4.4.7,  that  the  tuning  was  performed 
for  the  (7,6)  parameter  point.  The  criteria  used  for  tun¬ 
ing  was  to  choose  the  control  weighting  matrix  (U)  small 
(diagonal,  with  terms  of  3,  3,  and  9  along  the  diagonal) 
and  Increase  the  state  weighting  matrix  until  the  rms  error 


Table  5-5 


Compar i 
Contr 


son  of  Eigenvalues  of  State  Weighting  Matrix 
ol  Weighting  Matrix,  and  Cross  Term  Matrix 

State  Weighting  Matrix  X 
XI  -  0.98168 
X2  -  0.1615808 
X3  -  0.630921 
X4  -  0.001221 
X5  -  0.0012396 
X6  -  0.9891725 

Control  Weighting  Matrix  U 


U1 

-  1. 

70465 

U2 

-  0. 

0005257 

U3 

-  0. 

00063244 

Cross 

Term 

1  Matrix 

SI 

-  0. 

0 

S2 

-  0. 

0 

S3 

-  0. 

0 

S4 

-  0. 

0001252 

S5 

-  0. 

0000998 

on  the  true  states  stop  changing  drastically.  The  X  and  U 
matrices  were  tuned  for  the  (7,6)  parameter  point,  meaning 
that  the  filter  located  at  parameter  point  (7,6)  was 
matched  against  the  truth  model  located  at  parameter  point 
(7,6).  The  (7,6)  parameter  point  was  chosen  so  that  the 
parameter  point  would  be  near  the  center  of  the  space  but 
not  exactly  corresponding  to  a  fixed-banlc  parameter  point 
of  a  fixed-ban)<  algorithm  with  coarse  discretization:  the 
nine  filters  included  in  the  fixed-ban)<  were  located  at 
parameter  points  (1,1),  (1,5),  (1,9),  (5,1),  (5,5),  (5,9), 
(9,1),  (9,5),  and  (9,9).  The  final  values  of  the  X  and  U 
matrices  are  given  in  Table  5-6.  The  X  values  were  chosen 

Table  5-6 

State  Weighting  and  Control  Weighting  Matrices 
(for  the  (7,6)  parameter  point) 

(Units:  position  -  in;  velocity  -  in/sec) 


Position 
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Pos i t ion 
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(X2) 

-  7917.8 
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of 
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1 
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Velocity 

o  f 

Node 

2 

(X4) 

-  82.0 

Position 
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Hub 
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-  1. 

Velocity 

0  f 
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(X6  ) 
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1 

Actuator 

at 

Node 

1 

(Ul) 

-  3.0 

Actuator 

at 

Node 

2 

(U2) 

-  3.0 

Actuator 

at 

Hub 

(U3) 

-  9.0 

intentionally  large  as  compared  to  the  U  values,  thereby 
causing  a  large  emphasis  to  be  placed  on  quelling  the  state 
estimates  to  zero  while  not  strongly  caring  how  much 
control  energy  it  takes  to  accomplish  this.  The  impact  of 
tuning  the  X  and  U  matrices  will  be  addressed  further  in 
Section  5.9. 

5.9.  Controller  Study 

Once  the  X  and  U  values  were  determined,  a  controller 
study  was  performed  to  compare  the  performance  of  the 
moving-bank  controller  with  a  fixed-bank  controller  and  a 
benchmark  controller.  The  controller  used  in  this  study  is 
shown  in  Figure  5-25.  In  the  operation  of  this  adaptive 
algorithm,  the  parameter  estimate  locked  onto  the  true 
parameter  in  less  than  0.75  sec.  This,  coupled  with  the 
fact  that  a  dither  signal  was  being  applied  to  the  system 
for  0.5  sec,  was  the  basis  for  not  using  a  transition 
period  of  a  partially  adaptive  controller  as  described  in 
Section  2.4,  and  instead  going  straight  into  the  controller 
scheme  set  up  in  Figure  5-25.  Therefore,  the  remainder  of 
this  effort  was  performed  with  an  open  loop  dither  signal 
of  magnitude  300  and  frequency  30  rad/sec  being  applied  for 
the  first  0.5  sec.  The  magnitude,  frequency,  and  duration 
of  the  dither  signal  was  determined  by  trial  and  error. 

The  dither  signal  was  to  excite  the  system  but  not  to  act 
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as  the  major  forcing  function.  Following  the  period  of 
applying  the  open  loop  dither  signal,  multiple  model 
adaptive  control  was  applied  to  the  system.  The  complete 
details  on  the  steps  to  accomplish  this  study  are  given  In 
Section  4.4.8. 

In  this  study,  a  benchmark  controller  was  needed  to 
compare  with  the  moving-bank  and  fixed-bank  multiple  model 
adaptive  controllers.  This  benchmark  controller  was  estab¬ 
lished  by  using  a  single  filter  and  controller  located  at 
parameter  (7,6)  matched  against  a  truth  model  at  parameter 
(7,6).  The  controller  scheme  used  In  the  benchmark  con¬ 
troller  Is  given  In  Figure  4-1,  except  the  benchmark  con¬ 
troller  does  not  estimate  the  true  parameter  (a).  The 
benchmark  controller  performance  for  parameter  (7,6)  Is 
given  In  Figure  5-26.  This  figure  plots  the  mean  value 
(symbol  x),  the  mean  plus  1  sigma  (symbol  +),  and  the  mean 
minus  1  sigma  (symbol  y)  of  each  of  the  six  true  states  and 
three  control  Inputs  for  the  entire  length  of  the  runs. 

The  purpose  of  this  figure  Is  to  show  that  the  transients 
are  short  lived  and  stable.  Therefore,  in  the  remalner  of 
this  section  only  the  post-transient  performance  Is  given 
in  the  figures. 

In  order  to  compare  the  post-transient  performance  of 
the  controllers  effectively.  Figures  5-27  to  5-29  were 
constructed  exactly  like  Figure  5-26,  except  the  plots 
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started  at  t=  6.0  sec  in  order  to  remove  the  initial 
large-amplitude  transient  effects.  Also  recall  that  the 
center  of  the  moving  bank  is  always  started  at  the  para¬ 
meter  point  (5,5),  which  corresponds  to  the  center  of  the 
parameter  space.  Figure  5-27  contains  the  plots  for  the 
benchmark  controller.  Figure  5-28  contains  the  plots  for 
the  moving-bank  controller,  and  Figure  5-29  contains  the 
plots  for  the  fixed-bank  controller.  The  most  important 
states  to  compare  are  the  position  states  because  these 
states  give  an  indication  of  movement  of  the  truss  and  it 
is  this  movement  that  the  controller  is  to  quell  to  zero. 
The  figures  show  that  the  moving-bank  and  the  benchmark 
controllers  performed  nearly  identically  for  the  three 
position  states.  The  similarity  between  the  benchmark  con¬ 
troller  and  the  moving-bank  controller  was  expected  because 
the  moving-bank  algorithm  estimate  of  the  true  parameter 
locked  onto  the  actual  true  parameter  value.  After  about 
0.5  sec,  99.9  percent  of  the  probabilty  weighting  is  being 
applied  to  the  filter  that  corresponds  to  the  true 
parameter  value  (i.e.,  center  of  the  bank),  therefore  the 
moving-bank  algorithm  Is  performing  essentially  the  same  as 
the  benchmark  controller.  There  were  some  differences 
however  between  the  moving-bank  controller  and  the  bench¬ 
mark  controller.  The  benchmark  controller  performed  much 
better  for  the  state  corresponding  to  the  velocity  at 


node  2f  although  more  Importantly,  the  benchmark  controller 
required  much  less  control  energy  than  the  moving-bank 
controller.  This  is  probably  due  to  the  fact  that  the 
moving-bank  controller  does  have  to  go  through  some  tran¬ 
sient  effects  before  the  estimate  of  the  true  parameter 
locks  onto  the  true  parameter.  Note  also,  that  the 
velocity  states  are  not  quite  in  steady  state.  However, 
this  does  not  pose  a  problem  because  the  figures  clearly 
show  that  the  velocity  states  are  approaching  steady  state 
condition.  For  this  set  of  simulation  runs  the  fixed-bank 
controller  performed  (see  Figure  5-29)  much  worse  than  the 
benchmark  controller.  The  reason  for  this  poor  performance 
is  two-fold.  First,  it  was  initially  assumed  that  the 
fixed-bank  controller  would  put  most  of  its  probability 
weight  on  the  four  filters  that  surround  the  true  parameter 
point.  Instead,  it  put  nearly  all  of  its  probability 
weight  (l.e.  1)  on  one  filter  and  bounced  this  heavy 
weight  factor  nonsymmetr leal ly  throughout  the  nine  filters 
in  the  fixed  bank  configuration.  In  addition,  the  X  and  U 
matrices  were  tuned  for  the  (7,6)  parameter  point.  There¬ 
fore,  the  locations  in  parameter  space  where  the  fixed-bank 
controller  is  placing  its  heavy  weighting  is  not  tuned  for 
the  X  and  U  matrices. 

To  show  this  need  for  tuning  X  and  U  matrices  for 
every  parameter  point.  Figures  5-30  and  5-31  were 
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generated.  In  both  these  figures,  the  X  and  U  values 
obtained  for  the  (7,6)  parameter  point  (see  Table  5-4)  were 
used;  however,  the  true  parameter  point  was  (5,5).  As  can 
be  seen,  both  the  moving-bank  and  fixed-bank  controllers 
are  unstable.  Then  the  X  and  U  values  were  retuned  for  the 
(5,5)  parameter  point  (see  Table  5-7).  This  time  both  the 
moving-bank  and  fixed-bank  controllers  provided  a  stable 
system  (see  Figures  5-32  and  5-33).  Therefore,  the  X  and  U 
values  must  be  tuned  for  each  parameter  point.  This  mls- 
tunlng  of  certain  parameter  points  could  cause  the  poor 
performance  of  the  fixed-bank  controller  for  the  (7,6) 
parameter  point.  However,  the  coarseness  of  discretization 
in  the  fixed-bank  controller  Is  still  an  issue:  even  with 
separate  tuning  for  each  Implemented  parameter  point,  there 
is  still  no  single  elemental  filter/controller  that  is 
particularly  well  tuned  for  a  parameter  value  far  from  any 
single  coarsely  discretized  point  value. 

In  addition  to  showing  that  each  parameter  point 
needs  to  have  its  X  and  U  matrices  tuned.  Figure?  ‘ 

5-34  also  show  that  all  three  controllers  are  , 

nearly  identically.  This  is  due  to  the  fa  •  ••< 
true  parameter  happens  to  correspond  t  ^  • 
fixed  bank,  then  the  fixed-bank  alq  .  . 

true  parameter  locks  onto  the  a  *  a 
fact,  the  fixed-bank  alg  r  :  • ' n  •  . 


Table  5-7 

State  Weighting  and  Control  Weighting  Matrices 
(for  the  (5,5)  parameter  point) 

(Units:  Position  -  in;  Velocity  -  in/sec) 

Position  of  Node  1  (XI)  -  7€2,70.0 
Position  of  Node  2  (X2)  -  7917.8 
Velocity  of  Node  1  (X3)  -  508.5 

Velocity  of  Node  2  (X4)  -  82.0 

Position  of  Hub  (X5)  -  10.0 

Velocity  of  Hub  (X6)  -  1.0 

Actuator  at  Node  1  (Ul)  -  3.0 


Actuator  at  Node  2  {U2)  -  3.0 
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moving-bank  algorithm.  Therefore,  now  both  the  fixed-  and 
moving-bank  controllers  are  converging  to  the  benchmark 
controller  performance.  However,  when  the  true  parameter 
does  not  correspond  to  a  filter  in  the  fixed  bank,  then  the 
fixed-bank  algorithm's  estimate  of  the  true  parameter  jumps 
throughout  the  parameter  space  and  its  performace  never 
converges  to  the  performance  of  the  benchmark  controller, 
while  the  moving-bank  controller  characteristics  are 
converging  to  those  of  the  benchmark  controller  in  both 
cases . 

5.10.  Summary 

The  results  of  this  thesis  effort  were  presented  in 
this  chapter.  A  number  of  important  points  were  uncovered 
in  this  research.  First,  the  problem  with  the  moving-bank 
algorithm's  performance  in  last  year's  effort  by  Karnick 
was  due  to  the  values  of  the  measurement  noise  covariance 
and  dynamics  noise  strength.  Next,  the  space  discreti¬ 
zation  study  showed  that  the  parameter  space  discretization 
plays  a  large  role  in  the  performance  of  the  MMAE.  Follow¬ 
ing  the  space  discretization  study,  the  cross  term  study 
showed  that  adding  cross  terms  to  the  quadratic  cost 
function  did  not  substantially  improve  the  MMAC  perform¬ 


ance.  Finally,  the  moving-bank  MMAC  and  the  fixed-bank 
MMAC  were  compared  with  a  benchmark  controller.  The 
results  showed  that  the  performance  of  the  moving-bank 


MMAC  converged  to  that  of  the  benchmark  controller  because 
of  the  speed  at  which  the  estimated  true  parameter  con¬ 
verges  on  the  actual  true  parameter  value.  In  contrast, 
the  fixed-bank  algorithm  performance  converged  to  that  of 
the  benchmark  controller  if  the  true  parameter  happened  to 
correspond  to  a  filter  in  the  fixed  bank;  however,  if  the 
true  parameter  did  not  correspond  to  a  filter  in  the  fixed 
bank,  then  the  performance  was  substantially  worse. 
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Conclusions  and  Recommendations 


6.1.  Introduction 

The  Investigation  of  the  moving-bank  Multiple  Model 
Adaptive  Estimator  showed  very  encouraging  performance  in 
parameter  estimation,  state  estimation,  and  control  per¬ 
formance.  The  research  performed  in  the  effort  should 
provide  for  a  strong  base  for  follow-on  efforts. 


6.2.  Conclus ions 

The  values  of  the  measurement  noise  covariance  (R) 
play  an  extremely  large  role  in  the  performance  of  the 
moving-bank  algorithm.  In  fact,  this  was  the  predominant 
problem  with  the  moving-bank  algorithm  performance  In  last 
year’s  effort  191.  The  important  influence  that  R  had  on 
the  moving  bank  performance  was  expected  because  large  R 
values  would  tend  to  mask  the  differences  between  good  and 
bad  filters,  thereby  making  it  extremely  difficult  to  dis¬ 
cern  anything  distinctly  about  the  system  parameters.  Thl 
limit  on  the  R  values  is  the  major  driver  of  the  type  of 
measurement  devices  used  in  a  moving-bank  MMAE  problem. 

Once  the  moving  bank  is  being  pulled  in  the  right 
direction  through  proper  choice  of  accurate  enough  sensors 
additional  performance  can  be  achieved  by  determlng  the 
optimal  space  discretization  level.  The  correct 
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discretization  level  can  transform  a  poorly  operating 
moving  bank  into  an  effectively  operating  moving  bank. 

Again  this  is  not  surprising,  because  the  easier  it  is  for 
the  bank  to  distinguish  between  good  and  bad  filters,  the 
better  the  moving-bank  HMAE  operates,  both  in  state 
estimation  and  parameter  estimation. 

After  the  best  moving-bank  MHAE  has  been  designed, 
the  next  step  is  to  enhance  the  controller  performance  of 
the  MMAC.  For  the  large  spacestructure  model,  the 
moving-bank  MMAC  needs  to  have  the  state  weighting  matrix 
and  the  control  weighting  matrix  tuned  for  each  parameter 
point.  In  addition  to  the  tuning  of  the  weighting  matri¬ 
ces,  the  cross  terms  were  shown  to  have  little  effect  on 
the  quadratic  cost  function  for  this  problem's  choice  of 
sample  period.  Also,  the  fixed-bank  MMAC  performs  as  well 
as  the  moving-bank  MMAC  if  the  true  parameter  corresponds 
to  one  of  the  nine  filters  contained  in  the  fixed-bank 
MMAC,  but,  the  fixed-bank  MMAC  performance  degrades  sub¬ 
stantially  if  the  true  parameter  does  not  correspond  to  one 
of  the  filters  in  the  fixed-bank.  Thus,  the  moving-bank 
algorithm  inherently  has  greater  performance  potential, 
since  it  does  not  require  such  a  coarse  discretization  of 
parameter  space  in  order  to  yield  a  practically  imple- 
mentable  controller. 


6.3.  Recommendations 


It  is  strongly  recommended  that  research  continue 
using  the  rotating  two-bay  truss  model  of  a  flexible  space 
structure . 

The  most  important  recommendation  by  far  is  that  the 
effect  of  unmodelled  effects  be  Investigated.  This  would 
be  accomplished  by  keeping  the  filter  dimensioned  at  six 
but  increasing  the  "truth"  dimension  to  16  states  or  more. 
So  far  all  the  research  has  been  conducted  with  equal  order 
truth  and  filter  models.  This  is  not  very  realistic. 
Robustness  of  the  parameter  adaptation  process,  and  ulti¬ 
mately  of  the  state  estimation  and  control  per f ormanceSj 
requires  investigation. 

It  is  also  recommended  that  the  state  weighting 
matrix  and  control  weighting  matrix  be  tuned  for  each 
parameter  point.  This  would  entail  minor  modifications  to 
the  preprocessor  and  would  allow  the  moving-bank  MMAC  to 
obtain  better  control  performance. 

In  this  study  the  true  parameter  value  was  always 
held  constant,  however  the  parameters  can  also  be  slowly 
varying  or  make  jumps  in  parameter  space.  Therefore,  a 
study  needs  to  conducted  for  varying  true  parameter  values. 
The  performance  of  the  moving-bank  algorithm  could  be 
monitored  so  that  once  the  performance  degrades  to  a 
predetermined  pointy  then  a  reacqulsltlon  cycle  can  be 


performed  to  reacquire  a  good  estimate  of  the  true 
parameter . 

Finally,  a  performance  comparison  could  be  conducted 
between  both  the  state  estimates  and  the  resulting  controls 
produced  with  the  Bayesian  estimator  and  the  state 
estimates  and  controls  produced  with  a  MAP  estimator. 


Recall,  the  Bayesian  estimation  is  formed  by  a 
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Appendix  A:  LQG  Controller  Development  C9J 


Assume  the  following  stochastic  system  [3:20-22; 


6:33-35J 


x(t)  =  Fx(t)  +  Bu(t)  +  Gw(t) 


(A-1) 


where 


E[w(t)]  =  0 

E[^(t)y'^(t  +  r)]  =  Q  «(t) 

and  the  quadratic  cost  function  to  be  minimised  is: 


J  =  Ej/  (l/2)i:x'^(t)W  x(t)  +  u'^(t)w  u(t)3dt}  (A-2) 

I  'Nr 


where  w  and  w  are  weighting  matrices  to  be  chosen 
(iteratively)  to  yield  a  controller  with  desirable 
performance  characteristics.  The  optimal  constant-gain 
discrete  linear  feedback  control  law,  assuming  full-state 
access,  is  given  by: 


u(t^  )  =  -G^x(t^ 


(A-3) 


where  the  constant  gains,  G^,  that  minimize  J,  are  given 
by:  C23:68,1223 


(A-4) 


where  K 


and 


<*>  (t2,t 

and  the 
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satisfies  the  algebraic  Rlccatl  equation 
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is  the  state  transition  matrix  from  t^  to  t2 
<t>'s  in  Equations  (A-4)  and  (A-5)  are: 


<*>  =  4>  (t^,,,t.) 


It  should  be  noted  that  Equation  (A-3)  is  also  the 
solution  for  the  deterministic  LQ  optimal  control  problem 


with  no  driving  noise  w(t).  If  full  state  access  is 

replaced  by  noise-corrupted  measurements,  x(t^)  in 

Equation  (A-3)  is  replaced  by  the  state  estimate  ittt), 

which  is  often  described  ate  having  the  "certainty  equiva 
lence"  property  C21:171. 


Appendix  B:  Rotating  Two-Bay  Truss  System 
Matrices 


This  appendix  lists  the  system  matrices  for  the 
rotating  two-bay  truss.  The  reduced  order  matrices  are 
developed  from  the  mass  and  stiffness  matrices  {see 
Sections  3.2  and  3.4).  The  system  is  in  modal  formulation 
(see  Section  3.3)  and  is  composed  of  6  states  with  6 
measurements  and  4  control  Inputs.  The  six  states  were 
position  at  hub,  position  at  actuator  1,  position  at 
actuator  2,  velocity  at  hub,  velocity  at  actuator  1,  and 
velocity  at  actuator  2  respectively.  The  six  measurements 
were  position  at  actuator  1,  position  at  actuator  2, 
velocity  at  actuator  1,  velocity  at  actuator  2,  position  at 
hub,  and  velocity  at  hub  respectively. 


System  Matr ices 


St  1  £ f ness  Matrix 


ROW 

1 

1.188E+3 

1 . 966E+2 

.  OOOE  +  0 

.  OOOE  +  0 

! 

-6 . 424E  +  2 

.  OOOE  +  0 

-5. 461E+2 

-1 .966E+2 

.  OOOE+0 

.  OOOE+0 

. OOOE+0 

. OOOE+0 

ROW 

2 

1 .966E  +  2 

6 . 263E+2 

0 . OOOE+0 

-5.556E+2 

;• 

0  .  OOOE  +  0 

0 . OOOE+0 

-1.966E+2 

-7 . 077E+1 

> 

0 . OOOE+0 

0  .  OOOE  +  0 

O.OOOE+0 

0 . OOOE+0 

■ 

• 

Row 

3 

0 . OOOE+0 

0 . OOOE+0 

1 . 188E+3 

-1 .966E+2 

« 

M 

-5 . 461E+2 

1 . 966E  +  2 

-6. 424E+2 

0 . OOOE+0 

• 

0  .  OOOE  +  0 

0  .  OOOE  +  0 

0 . OOOE+C 

0 . OOOE+0 

m 

Row 

4 

0  .  OOOE  +  0 

-5, 556E+2 

-1.966E+2 

6 . 263E+2 

• 

i 

1 . 966E  +  2 

-7 . 077E+1 

0 . OOOE+0 

0 . OOOE+0 

0 . OOOE+0 

0 . OOOE+0 

O.OOOE+0 

0  .OOOE  +  0 

Row 

5 

-6, 424E+2 

0 . OOOE+0 

-5.461E+2 

1.966E+2 

U 

4 . 019E+3 

6 .69  3E  +  1 

0 . OOOE+0 

0 . OOOE+0 

-2.099E+3 

0 . OOOE+0 

-7.320E+2 

-2.635E+2 

■ 

Row 

6 

O.OOOE+0 

0 . OOOE+0 

1 .966E+2 

-7.077E+1 

6.693E+1 

7.212E+2 

O.OOOE+0 

-5. 556E+2 

N 

0 . OOOE+0 

0 . OOOE+0 

-2.635E+2 

-9  .  487E  +  1 

* 

Row 

7 

-5. 461E+2 

-1.966E+2 

-6.424E+2 

0  .  OOOE  +  0 

■ 

0  .  OOOE  +  0 

0 . OOOE+0 

4 . 019E+3 

-6 .693E+1 

y 

-7. 320E+2 

2 .635E  +  2 

-2.099E+3 

0 . OOOE+0 

« 

y 

Row 

8 

-1 .966E  +  2 

-7 . 077E  +  1 

0 . OOOE+0 

0  .  OOOE  +  0 

y 

0 . OOOE+0 

-5.556E+2 

-6.693E+1 

7 . 212E+2 

\ 

2 . 635E  +  2 

-9 . 487E+1 

O.OOOE+0 

0  .  OOOE  +  0 

• 

Row 

9 

0 . OOOE+0 

0 . OOOE+0 

O.OOOE+0 

0  .  OOOE  +  0 

-2 .099E+3 

0 , OOOE+0 

-7.320E+2 

2 .635E+2 

8 . 618E  +  4 

4 . 788E+4 

0. OOOE+0 

0 . OOOE+0 

Row 

10 

0 . OOOE+0 

0 . OOOE+0 

O.OOOE+0 

0 . OOOE+0 

1 

0 . OOOE+0 

0 . OOOE+0 

2.635E+2 

-9 . 487E+1 

4.788E+4 

1 . 390E+5 

O.OOOE+0 

-1 . lllE+5 

Row 

11 

O.OOOE+0 

0 . OOOE+0 

O.OOOE+0 

0  .  OOOE+0 

-7. 320E+2 

-2, 635E+2 

-2.099E+3 

0 . OOOE+0 

■s 

0 . OOOE+0 

0 . OOOE+0 

8.618E+4 

-4 . 788E+4 

i 

c 

Row 

12 

0 . OOOE+0 

0 . OOOE+O 

0 . OOOE+0 

0 . OOOE+0 

• 

-2 . 635E  +  2 

-9 . 487E+1 

O.OOOE+0 

0 . OOOE+0 

% 

V 

0 . OOOE+0 

-1 . lllE+5 

-4 .788E+4 

1 . 390E  +  5 

\ 

Mass  Matrix 

Note  that  the  first  8  diagonal  elements  are 
essentially  the  values  of  the  non-structural  mass  because 
the  non-structural  mass  is  large  compared  to  the  structural 
mass  . 
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0 . OOOE+0 
-8 , 152E-7 

0  .  OOOE+0 

0 . OOOE+0 

1 . 294E+0 

0 . OOOE+0 

-2.395E-6 

O.OOOE+0 

3.210E-6 

8 . 621E-7 
7 .764E-7 
1.156E-6 

Row 

7 

6 . 652E-6 

0  .  OOOE  +  0 
8. 916E-6 

2.395E-6 

0 . OOOE+0 
-3.210E-6 

6.927E-6 

1.294E+0 

2.263E-5 

0 . OOOE+0 
8.152E-7 
0  .  OOOE  +  0 

Row 

8 

2 . 395E-6 

0 . OOOE+0 
-3 . 210E-6 

8 . 621E-7 
7.764E-7 

1 . 156E-6 

0 . OOOE+0 
8.152E-7 
O.OOOE+0 

0  .  OOOE  +  0 
1 . 294E  +  0 
0  .  OOOE  +  0 

Row 

9 

0 . OOOE+0 
2.263E-5 

8 . 817E-4 

0 . OOOE+0 

0 . OOOE+0 
-6 . 402E-5 

0  .  OOOe  +  0 

8 .916E-6 

0 . OOOE+0 

0 . OOOE+0 
-3.210E-6 
0 . OOOE+0 

Row 

10 

0 . OOOE+0 

0 . OOOE+0 
-6 .402E-5 

0 . OOOE+0 

0 . OOOE+0 

8 . 343E-4 

O.OOOE+0 

-3.210E-6 

O.OOOE+0 

0  .  OOOE  +  0 
1.156E-6 
1 . 553E-4 

Row  11 


0 . OOOE+0 
8.916E-6 
0 . OOOE+0 


0 . OOOE+0 
3. 210E-6 
0 . OOOE+0 


0 . OOOE+0 
2.263E-5 
8 . 817E-4 


0 . OOOE+0 
0 . OOOE+0 
6 . 402E-5 


t "  I 

•C*; 

I 


Row  12 

0 . OOOE+O 

0 .OOOE+O 

0. OOOE+O 

0 .  OOOE+O 

3.210E-6 

1 . 156E-6 

0. OOOE+O 

0  .  OOOE  +  O 

i 

0. OOOE+O 

1.553E-4 

6.402E-5 

8 . 343E-4 

VA 

A  Matrix 

•rf; 

0 , OOOE+0 

0 .OOOE+O 

0  .OOOE+O 

1. OOOE+O 

0. OOOE+O  0. OOOE+O 

O.OOOE+0 

0, OOOE+O 

0. OOOE+O 

0, OOOE+O 

1. OOOE  +  O  0. OOOE  +  O 

0 . OOOE+0 

0 .OOOE+O 

0 .OOOE+O 

0. OOOE+O 

0. OOOE  +  O  1. OOOE  +  O 

0 . OOOE+0 

0 . OOOE+O 

0  .  OOOE  +  O 

0 . OOOE+O 

0, OOOE+O  0. OOOE  +  O 

0. OOOE+0 

-7 ,918E  +  1 

0 .OOOE+O 

0. OOOE+O 

-8.898E-2  0. OOOE+O 

0 . OOOE+O 

0  .  OOOE  +  O 

-5 . 085E  +  2 

0 . OOOE+O 

0. OOOE+O  -2.255E-1 

1 

Rows  1  and  4 

correspond 

to  the  rigid  body  angular 

/v 

position  and  velocity,  respectively.  Rows  2  and  3 

•v 

Cv! 

correspond  to  the  position  of  the  first  and  second  bending 


modes;  respectively,  while  rows  5  and  6  represent  the 
velocity  of  these  bending  modes,  respectively. 

« 

I  .• 

B  Matrix 


0 . OOOE+0 
0 . OOOE+0 
0 . OOOE+0 
-4 .773E-1 
2 . 589E-1 
9 . 488E-2 


0 . OOOE+0 
0 . OOOE+0 
0 . OOOE+0 
-2.711E-1 
-4.736E-1 
1 . 590E-1 


O.OOOE+0 
0 . OOOE+0 
O.OOOE+0 
1 . OOOE+0 
O.OOOE+0 
O.OOOE+0 


1 . OOOE+0 
0 . OOOE+0 
0 . OOOE+0 
0 . OOOE+0 
0 . OOOE+0 
0 . OOOE+0 


The  first  two  columns  represent  inputs  from 
actuators  located  on  the  truss  while  the  third  and  fourth 
columns  are  due  to  an  actuator  located  on  the  hub  which 
induces  inputs  at  both  the  position  and  velocity  level  (see 


Section  4.4.4.),  The  non-zero  portion  of  the  first  two 
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columns  was  designed  in  physical  coordinates  In  the 
unreduced  system  as 

T 

01000000000  o] 

P  0000100000  oj 

and  then  transformed  Into  modal  coordinates.  The  angular 
input  actuator  was  then  added  by  augmenting  a  column  and 
entering  a  1.00  In  the  row  corresponding  to  the  angular 
velocity  state  and  the  angular  position  state. 

H  Matrix 

6.066E-1  -3.239E-1  -1.228E-1  O.OOOE+0  O.OOOE+0  O.OOOE+0 

3.444E-1  5.925E-1  -2.057E-1  O.OOOE+O  O.OOOE+0  O.OOOE+0 

O.OOOE+O  O.OOOE+O  O.OOOE+O  6.066E-1  -3.239E-1  -1.228E-1 

O.OOOE+O  O.OOOE+O  O.OOOE+O  3.444E-1  5.925E-1  -2.057E-1 

l.OOOE+0  O.OOOE+O  O.OOOE+O  O.OOOE+O  O.OOOE+O  O.OOOE+O 

|0.000E+0  O.OOOE+O  O.OOOE+O  l.OOOE+0  O.OOOE+O  O.OOOE+0_ 

The  first  and  second  rows  represent  measurements 
from  the  position  sensors  located  at  node  1  and  node  2. 

The  third  and  fourth  rows  represent  measurements  from 
velocity  sensors  located  at  node  1  and  node  2.  Finally, 
the  last  two  rows  represent  measurements  of  the  angular 
position  and  velocity  of  the  hub  (see  Section  3.4.4.).  The 
position  and  velocity  portion  of  the  matrix  were  calculated 
in  physical  coordinates  and  then  transformed  into  modal 
coordinates  (see  Section  3.3.),  where  the  angular 


measurements  were  augmented.  Note  that  the  entries  in  rows 
1  and  2  are  identical  to  those  in  rows  3  and  4  because  of 
co-location  of  position  and  velocity  sensors.  Rows  1  and  2 
were  designed  in  physical  coordinates  in  the  unreduced 
system  as : 


010000000000 
0  00001000000 


Only  the  first  twelve  columns  are  listed  because  the 
remaining  columns  are  zero. 


D 


Matrix 


-4 . 253E-4 
6 . 396E-5 
O.OOOE+0 
0 . OOOE+0 
0 . OOOE+0 
0 . OOOE+O 


6 . 396E-5 
-3 . 583E-4 
0 . OOOE+O 
0 . OOOE+O 
0 .OOOE+O 
0 .OOOE+O 


0 . OOOE-0 
0 . OOOE-0 
0 . OOOE+O 
0. OOOE+O 
0 .OOOE+O 
0 . OOOE+O 


0 .OOOE+O 
0 . OOOE+O 
0 . OOOE+O 
0 . OOOE+O 
0 . OOOE+O 
0 . OOOE+O 


The  matrix  (see  Section  3.5)  is  used  in  a 
method  of  order  reduction  referred  to  as  singular 
perturbations.  This  thesis  used  reduced  order  models; 
however,  this  method  was  not  implemented  because  truth 


models  were  of  the  same  order. 
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Abstract 

This  investigation  focused  on  the  use  of  moving-bank 
multiple  model  adaptive  estimation  and  control  (MMAE). 
Moving-bank  MMAE  reduces  the  computational  burden  of  MMAE 
by  implementing  only  a  subset  of  the  Kalman  filters  (.^ 
filters  versus  100  in  this  research)  that  are  necessary  to 
describe  the  system  to  be  estimated/controlled.  Important 
to  the  development  of  the  moving-bank  MMAE  are  the  derision 
logics  governing  the  selection  of  the  subset  of  filters. 

The  decision  logics  cover  three  situations:  initial 
acquisition  of  unknown  parameter  values;  tracking  unknown 
parameter  values;  and  reacqulsltlon  of  the  unknown  prar- 
meters  following  a  "jump"  change  in  these  parameter  values. 

This  thesis  applied  moving-bank  MMAE  to  a  rotating  two 
bay  truss  model  of  a  flexible  spacestructure .  The  rotating 
two  bay  truss  approximated  a  space  structure  that  had  a  hub 
with  appendages  extending  from  the  structure.  The  mass  of 
the  hub  was  large  relative  to  the  mass  of  the  appendage. 

The  mathematical  model  was  developed  using  finite  element 
analysis,  transformed  into  i.iOdal  formulation,  and  reduced 
using  a  method  referred  to  as  singular  perturbations. 
Multiple  models  were  developed  by  assuming  that  variation 
occurred  in  the  mass  and  stiffness  of  the  structure. 

Results  showed  that  the  primary  cause  for  the  Increase 
in  the  performance  of  the  moving-bank  algorithm  over  past 
efforts  was  the  choice  of  the  measurement  covariance  noise 
matrix.  This  research  effort  also  showed  that  the  space 
discretization  level  can  play  an  Important  role  in  the 
performance  of  the  moving-bank  MMAC;  therefore,  care  must 
be  taken  in  determining  the  discretization  level.  Finally, 
this  research  effort  showed  that  the  moving-bank  MMAC 
provided  increased  control  performance  over  a  fixed  bank 
controller  with  a  discretization  that  covered  the  entire 
range  of  parameter  variation.  In  addition,  the  moving-bank 
controller  converged  to  the  performance  of  a  benchmark 
controller . 
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